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Preface  and  contents 

The  three  main  directions  of  the  research  work  conducted 
by  the  VLP  team  at  Innsbruck  daring  the  past  four  years 
include : 

1.  Investigation  of  VLP  propagation  through  solid  media, 
preferably  rock,  and  construction  of  the  necessary 
measuring  apparatus . 

2.  Measurement  of  the  electrical  conductivity  S' and  the 
dielectric  constant  £  of  various  rock,  preferably  in 
natural  state. 

3.  Theoretical  studies  on  problems  of  excitation  and 
propagation  of  VLP  waves  in  solid  media. 

In  the  present  report,  the  research  results  of  1964  are 
described.  For  a  better  survey  and  especially  for  better 
orientation  of  laymen  in  the  field  of  VLP  propagation,  a 
brief  review  shall  follow. 

On  the  first  main  direction  (chapters  1  and  2): 

In  chapter  1  (Apparatus)  a  3mall,  easily  portable  field 
strength  meter  for  frequencies  ranging  from  1.9  to  25  kc/sec 
is  described.  The  nuvistorized  apparatus  has  a  high  sensiti¬ 
vity  and  a  large  signal-to-noise  ratio.  A  wire-wound  ferrite 
rod,  approximately  1  m  long  is  used  as  an  antenna. 

The  studies  on  optimum  receiving  antennas  for  VLP  signals 
have  been  terminated.  The  dimensioning  of  magnetic  dipole 
antennas  with  ferrite  cores  is  discussed  in  detail. 

In  the  mine  of  St.  Gertraudi,  a  new  transmitting  antenna 
with  an  antenna  area  of  l600  m2  and  a.  wire  cross  section  of 
135  mm2  was  set  up  in  the  immediate  neighborhood  of  the  VLP 
laboratory  there.  The  distances  reached  with  this  new  antenna 
are  much  larger  than  those  reached  with  former  antennas. 

For  studying  the  exact  direction  of  the  electromagnetic 
field  of  a  VLP  transmitter,  a  direction-finder  antenna  was 
constructed  that  can  be  freely  rotated  in  space  and  fixed 
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in  any  desirable  position.  The  results  of  measurement  are  dis¬ 
cussed  in  chapter  1. 

For  exact  angular  measurements  in  other  mines  when  using 
the  portable  transmitter,  a  small  transmitting  antenna  w as 
designed.  Any  angular  adjustment  becomes  possible  with  this  an¬ 
tenna,  up  to  an  accuracy  of  +2°.  It  is  mainly  used  for  deter¬ 
mining  tho  directional  patterns  and  field  directions  in  the 
neighborhood  of  geological  interference  zones. 

In  chapter  2,  the  propagation  measurements  in  homogeneous 
and  in  inhomogeneous  rock  of  St.  Gertraudi  and  Bleiberg,  re¬ 
spectively,  are  dealt  with.  A  50  SI  transmitter  was  used  for 
measuring  the  propagation  and  the  directional  pattern.  The 
attenuation  values  thus  obtained  cannot  be  compared  with  the 
attenuation  of  a  remote  VLF  transmitter  (GBR).  It  is  shown 
that  the  great  differences  in  attenuation  are  assumed 
be  due  to  the  fact  that  one  measurement  was  conducted  in  the 
near  field;  whereas  the  other  one  was  made  in  the  far  field. 

The  difference  in  the  type  of  propagation  and  attenuation  ob¬ 
served  in  the  above  two  cases  can  also  be  derived  theoretically. 

The  probler  of  near  field  and  far  field  is  discussed  in 
another  section  giving  a  number  of  diagrams  of  the  field  strength 
trend  as  dependent  on  the  quantities  r,^,  u>  and  6". 

Measurements  in  inhomogeneous  media  made  it  necessary  to 
determine  the  value  of  the  field  strength  at  the  point  of 
reception  as  well  as  the  direction  of  the  field  vector.  If  the 
rock  contains  inhomogeneities,  the  latter  may  be  determined 
much  more  easily  and  accurately  by  measuring  the  field  direction 
than  by  measuring  only  the  field  strength  value. 

On  the  second  main  directio.  (chapter  J>)i 

Measurements  with  direct  current  and  with  alternating 
current  wore  conducted  the  latter  in  samples  as  well  as  in 
solid  rock. 

By  means  of  measurements  following  the  yenner  method  it 
was  empirically  shown  from  what  electrode  distance  onward  the 
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effect  of  the  gallery  space  on  the  result  of  measurement  may 
be  neglected.  It  can  also  be  seen  that  the  Wenner  method  is  a 
four-electrode  configuration  well  suited  for  measurements  in 
mine  galleries  (section  3*l)« 

The  sample  measurements  yielding  a  strong  frequency  de¬ 
pendence  of  £,  and  6  were  evaluated  by  K.W.  Wagner's  dispersion 
theory,  assuming  a  continuous  distribution  of  relaxation  times. 
The  graphic  methfd  by  W.A.  Yager  was  used  for  the  evaluation, 
although  neither  the  results  obtained  for  completely  dry  nor 
those  obtained  for  completely  damp  samples  can  be  evaluated  by 
means  of  the  curves  calculated  in  his  studies.  The  evaluable 
results,  however,  show  a  good  agreement  between  experiment 
and  theory.  The  dependence  of  the  humidity  constants  contained 
in  the  theory  was  studied  (section  3*2). 

Furthermore  the  propagation  resistance  of  electrodes  in 
solid  rock  ns  dependent  on  the  frequency  was  measured.  The 
Rvalues  calculated  therefrom,  however,  could  not  be  made  con¬ 
sistent  with  the  results  obtained  for  the  samples  (section  3»3)» 

On  the  third  main  direction  (chapter  4): 

The  theoretical  part  of  the  present  annual  report  is  mainly 
dedicated  to  one  region  essential  for  the  excitation  of  VLF 
waves;  Dipole  radiation  in  a  homogeneous  and  conducting  medium. 
So  far,  dipole  radiation  in  air  or  in  vacuo  as  an  elementary 
process  of  every  antenna  phenomenon  has  been  treated  adequately 
in  many  publications,  and  has  been  investigated  in  great 
detail.  Furthermore  it  has  been  pointed  out  that  the  mathe¬ 
matical  representation  that  describes  the  physical  process 
when  passing  from  vacuum  to  a  conducting  medium  must  be  inter¬ 
preted  with  great  care,  and  that  for  certain  physical  standard 
concepts,  the  definitions  used  for  processes  in  vacuo  have 
to  be  modified  for  describing  the  processes  in  a  conducting 
medium.  So  far,  the  relation  between  the  two  types  of  dipoles, 
electric  and  magnetic,  has  not  been  studied  sufficiently  as 
regards  their  energy  balance.  Therefore,  a  distinct  comparison 
between  the  electric  and  the  magnetic  dipoles  for  VLF  waves  in 
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conducting  media  has  been  described  in  section  (4*2). 

Besides  the  above  problems,  it  was  our  task  to  develop  the 
theory  of  conductivity  measurements  by  means  of  the  potential 
distribution  of  steady  currents  and  thus  give  an  estimation 
to  what  an  extent  the  experiments  conducted  in  mi.te  galleries 
below  ground  can  be  described  by  the  usual  expressions  that 
hold  for  homogeneous  full  spaces  or  semispaces. 
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1 .  Apparatus 

1.1  Nuvistorized  field  strength  meter  E4  for  the  VLF  range 

Introduction 


A  number  of  field  strength  meters  were  designed  within  the 
framework  of  the  VLF  research  program,  and  experiments  with  various 
antenna  designs  were  conducted  (see  annual  reports).  Owing  to  the 
range  of  application  -  measurements  in  mines  -  the  common  character¬ 
istics  of  the  designs  were  small  dimensions  of  the  measuring  device 
as  well  as  the  antenna,  and  operation  independent  of  the  mains  supply. 

When  developing  the  field  strength  meter  E4,  new  means  of  increasing 
the  sensitivity  and  signal-to-noise  ratio  were  found. 


Design  characteristics  of  the  field  strength  meter  E4 

The  antenna  of  a  field  strength  meter  for  very  low  frequencies 
will  in  any  case  be  small  (its  electric  length  as  well  as  its 
geometrical  dimensions)  as  compared  to  the  wave  length  A.  .  In  the 
measuring  range  of  the  described  device,  A.  is  about  10  to  500  km 
in  air.  Coupling  of  a  ’’short"  antenna  to  the  transmitter  field  will 
be  correspondingly  loose.  The  coupling  factor  will  be  increased  on 
the  one  hand  by  increasing  .he  area  of  a  frame  antenna,  and  on  the 
other  hand  by  using  highly  permeable  ferrite  cores  by  means  of  which 
the  lines  of  force  in  the  antenna  coil  can  be  concentrated.  The  use  of 
large  frame  coils  as  antennas  had  to  be  abandoned  after  several 
experiments,  as  the  attainable  Q  values  were  far  too  low.  In  order 
to  reach  an  induced  voltage  U  as  high  as  possible  when  using  inductive 
antennas  with  ferrite  cores  at  a  given  field  strength  E,  the  effective 
height  hg££  which  is  very  small  for  this  antenna  type,  has  to  be 
increased  (see  the  relation  U  =  he£f  •  E  ).  It  must  therefore  be 
the  aim  of  the  investigation  and  development  to  obtain  a  low-loss 
antenna  circuit  which  has  to  be  coupled  to  an  amplifier  with  high 

-i  V\  rui  4*  1  mno/l  OV%/>  A  O  «-»  A  4  4-  J  J  1  AMMAM  -I  M  i.  A  t.  . 
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kept  small  as  compared  to  the  radiation  damping  caused  by  the  radiation 
resistance  oT  the  antenna. 
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These  requirements,  however,  can  hardly  be  fulfilled  with 
transistorized  amplifiers,  since  their  comparatively  low-resistance 
input  damps  the  antenna  circuit  even  at  very  accurate  matching. 
Furthermore  it  is  impossible  to  utilise  the  v/hole  circuit  energy 
for  controlling  the  transistor,  especially  when  using  a  series- 
resonant  circuit. 

A  number  of  studies  were  also  dedicated  to  the  noise  of  the 
input  stage,  because  the  signal-to-noise  ratio  should  be  as  high 
as  possible  even  at  small  received  field  strengths. 

The  use  of  highly  permeable  ferrites  as  core  material  for  the 
antenna  rods  makes  possible  the  construction  of  parallel  resonant 
circuits  with  a  high  Q  value.  High  Q  values  of  the  circuit  decrease 
the  band  width  and  increase  the  signal-to-noise  ratio  k,  which 
is  gi/en  by 
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defined  for  an  electrical  field  strength  E  =  1. 

ug  =  signal  voltage  ^eff  =  e^’^ec'tlve  antenna  height 

Ur  =>  noise  voltage  n  =  number  of  turns  of  the  anten- 

Q  =  quality  of  circuit  na  coil 

p  =  permeability  at  the  coil  ^  =  received  wavelength 

Or 

location  having  a  distance  a  from  the  rod  center 
=  B  /h  (B  =  magnetic  induction  at  the  point  of  location) 

=  Pef£  =  effective  permeability  when  the  coil  is  placed  in  the 
rod  center. 

The  maximum  input  voltage  that  can  be  reached  at  a  given  field  strength 
is  as  follows; 

27i*n*pe££*E»F 


(1,2) 


max 


F  =  cross-sectional  area  of  the  coil. 

The  equations  1,1  and  1,2  show  that  it  is  d^sirabla  to  reach  high  n, 

Pa  ^efP*  an*  vplues* 

Details  on  the  use  of  ferrite  antennas  in  the  VLF  range  and 
their  calibration  shall  follow  in  a  separate  report. 

The  transistors  were  replaced  by  nuvistors  in  order  to  make  full 
use  of  the  high  Q  values  of  the  antenna  circuit  reached  by  means  of 
a  suitable  design,  i.e.  to  match  the  amplifier  to  the  high  resonant 
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resistance.  The  first  amplifier  stage  is  a  cascode  with  negative  feedback. 
The  equivalent  noise  resistance 

T 

Re  =  0.64  *  Y  °  gTs  ^  2‘5/S  ^  400  Q 

is  determined  only  by  the  nuvistor  triode  whose  cathode  is  grounded* 

It  is  negligibly  small  as  compared  to  the  resonant  resistance  Rq  which 

r~  . . — ■  1  ~ 

determines  the  noise  voltage  across  the  input:  =  V4KT^f»Ro  •  The 

cascode  stage  and  the  tuning  capacitors  were  directly  attached  to  the 
antenna  rod  in  order  to  reduce  screening  as  far  as  possible  and  avoid 
additional  circuit  damping.  A  multiple  cable  connects  this  antenna 
amplifier  with  the  second  part  of  the  measuring  device,  which  comprises 
a  selective  amplifier,  an  indicator,  and  current  supply. 

The  second  selective  circuit  was  no  longer  built  with  RC  networks, 
but  v/ith  LC  series  resonant  circuits.  High  Q  values  of  the  circuit  can 
be  reached  by  connecting  these  circuits  between  two  cathode  followers. 

The  utilizable  resonance  rise  of  the  signal  voltage  makes  another 
amplifier  stage  superfluous. 

Description  of  the  circuit  diagram  (Fig.  1,1) 

We  have  already  mentioned  that  the  field  strength  meter  E4 
consists  of  two  parts: 

I.  The  antenna  part  (left  of  the  dashed  line  in  the  circuit  diagram) » 
The  ferrite  rod  antenna  A  has  a  number  of  cross  coil  windings 
which  can  be  interconnected  by  means  of  the  switch  yielding  the 
following  values  of  inductance:  in  position  1  ...  2.19  H 

2  ...  600  mH 

3  ...  352  mH 

Resonant  frequencies  ranging  from  1.9  to  25  kc/sec  can  be  adjusted  by 
means  of  switch  S^,  the  nonvariable  capacitors  C^,  C^,  and  the 
variable  capacitors  C^,  C^,  Cg.  Range  overlapping  is  considerably  high 
owing  to  the  high  L-C  ratio.  The  socket  Bu^  is  used  for  feeding  a 
calibration  voltage,  and  for  further  experiments. 

The  alternating  voltage  induced  in  the  antenna  circuit  reaches  the 
input  of  a  cascode  amplifier  that  consists  of  the  nuvistors  Ro.  1  and 
Ro.2,  In  order  to  reach  a  maximum  resonance  rise  of  the  antenna  circuit, 
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the  latter  is  to  be  damped  but  slightly  by  the  tube  input.  Measurements 
yielded  an  optimum  grid  bias  of  1.35  V  which  corresponds  to  a  maximum 
tube  input  resistance.  The  grid  bias  is  produced  as  voltage  drop  at  the 
resistors  and  R^  by  the  common  anode  current.  The  resistor  R^  causes 
negative  feedback*  The  anode  circuit  of  the  cascode  is  coupled  with  a 
cathode  follower  (Rb*3)  acting  as  an  impedance  transformer  via  the 
high  pass  C^q  -  R^  to  attenuate  interfering  technical  frequencies. 

The  antenna  unit  is  pivoted  on  a  stand  and  connected  with  the 

II.  measuring  part 

via  a  screened  multiwire  cable  K.  In  order  to  keep  the  band  width  of  the 
measuring  device  small,  the  signal  is  fed  to  another  selective  circuit 
which  is  constructed  as  a  series  resonant  circuit.  A  range  switch 
consisting  of  4  switch  segments  serves  for  switching  over  on  the 
one  hand  the  inductances  through  L^,  and  on  the  other  hand  the 
capacitors  C^,.  through  V  In  one  position,  the  selective  circuit 
is  switched  off  whereby  the  band  width  of  the  unit  is  increased.  The 
cathode  follower  (R0.4)  connected  after  the  selective  circuit  provides 
low  circuit  damping  as  well  as  low  output  impedance.  The  adjusting 
potentiometer  is  used  for  amplification  calibration,  the  socket  Bu^ 
serves  as  output  (connection  of  an  oscilloscope  etc.).  Amplification 
may  be  attenuated  in  stages  (3  :  l)  by  means  of  a  calibrated  voltage 
divider  and  the  switch  S^.  The  second  amplifier  of  the  measuring  device 
is  the  nuvistor  stage  Ro.5.  The  amplified  output  signal  is  fed  a 
voltage  multiplier  stage  (diodes  1^  -  D^).  A  direct  voltage  drop 
proportional  to  the  signal  voltage  occurs  across  the  resistor  Rgg* 
direct  voltage  controls  the  anode  current  and  thus  also  the  length 
of  the  luminescent  band  of  the  tuning- indicator  tube  R0.6.  This  tube  has 
a  double  task:  it  makes  the  identification  of  keyed  transmitters  and 
strong  interferences  much  easier,  and  the  instrument  J  calibrated  in 
millivolt  input  voltage  lies  in  its  anode  circuit.  The  nonlinear  division 
of  the  scale  of  this  instrument  which  facilitates  the  reading  of  small 
voltages,  and  which  corresponds  to  the  characterista cs  of  the  tube,  is 
highly  advantageous.  If  there  is  no  input  signal,  the  instrument  pointer 
is  on  the  right-hand  side  of  the  scale,  corresponding  to  the  maximum 
tube  current,  i.e.  at  the  aero  point  or  the  beginning  of  the  mv  scale 
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(adjustment  by  the  trimming  potentiometer  with  the  switch 
being  in  its  lowest  position).  Overdriving  of  the  instrument  by 
extremely  strong  signals  i3  impossible,  since  the  instrument 
current  will  then  become  zero.  Furthermore  it  is  possible  to  use  the 
instrument  for  measuring  the  filament  voltage  and  the  anode  voltage 
by  means  of  the  switch  S^.  The  filament  voltage  is  indicated  on  an 
expanded  scale,  i.e„  the  zero  point  of  the  instrument  is  shifted  to 
4  v  by  means  of  the  Zener  diode  .  Thus,  a  slight  drop  in  filament 
voltage  can  be  read  much  more  easily;  this  is  of  special  importance  for 
accurate  field  strength  measurements.  Another  improvement  of  the 
stability  is  that  the  filament  voltage  of  the  indicator  tube  is 
kept  constant  by  means  of  a  stabilization  circuit  (transistors  Tr.  1 
and  Tr.  2,  Zener  diode  Dg). 

A  lead  storage  cell  serves  as  a  current  source  for  filament 
voltages,  and  an  NC  battery  for  the  anode  voltages.  The  switches 
Sg  and  are  used  as  operating  switches  for  measuring  and  charging.  An 
easily  portable  case  with  the  dimensions  230~x  130  x  140  mm  and  a  weight 
of  5*2  kg  contains  the  selective  amplifier,  indicator  and  current 
source  including  the  battery.  The  small  dimensions  of  the  device 
also  make  possible  measurements  in  difficult  terrain. 

Measured  values 

The  resonance  curve  of  the  antenna  unit  -  antenna  circuit  + 
nuvistor  cascode  +  cathode  follower,  without  the  selective 
amplifier,  is  shown  in  Fig,  1,2.  The  solid  curve  holds  for  a 
resonant  frequency  fQ  of  10  kc/sec,  the  dashed  line  holds  for 
3  kc/sec.  The  measuring  arrangement  for  taking  the  curves  is 
shown  in  the  block  diagram.  The  measuring  frequency  f  is  produced 
by  an  RC  generator  and  measured  with  a  digital  frequency  meter.  A  frame 
coil,  30  cm  in  diameter,  was  used  as  transmitting  antenna.  The  antenna 
stage  of  the  field  strength  meter  was  sot  up  at  a  distance  of  6  meters, 
a  precise  tube  voltmeter  being  connected  to  it. 

Measured  band  width  at  3  dB  voltage  drop: 

b  =4f  =  d  •  fQ  ...  at  5  kc/sec  :  ]0  cps 

at  10  kc/sec  :  75  cps 
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This  corresponds  to  a  value  Q  =  l/d  =  300  at  3  kc/sec 

=  133  at  30  kc/ sec 

The  hand  width  and  resonance  curve  of  the  field  strength  meter  with  a 
selective  amplifier  at  10  kc/sec  corresponds  approximately  to 
the  3  kc/sec  curve  of  Fig.  1,2. 

Selective  amplification  V  =  Uo/Le  at  10  kc/sec  :  1420-fold 

at  3  kc/sec  :  1580-fold 

(Ue  connected  to  Bu^,  Uo  is  measured  at  Bu2>  see  Fig.l) 

Smallest  measurable  signal  voltage  Ue  (at  Bu^)  :  5  pv  (accuracy  of 

instrument  reading) 

Maximum  measurable  signal  voltage  with  connected  selective 

amplifier  :  .  0  mv 
without  selective  amplifier  : 15C  mv. 

The  background  noise  of  the  field  strength  meter  with  antenna 
could  not  be  measured  as  there  is  no  suitable  Faraday  screen  for 
low  frequencies  available.  A  measurement  in  the  Schwaz  mine  was  con¬ 
ducted  to  get  an  idea  of  the  height  of  the  noise  level.  The  site  of 
measurement  b  ;low  ground  was  2.2  horizontal  km  away  from  the 
gallery  entrance.  Tho  maximum  depth  was  1000  m.  The  noise  signal 
occurring  at  the  output  Bu2  was  measured  by  an  oscilloscope. 

This  noise  is  composed  of  the  background  noise,  the  local  noise  level, 
and  atmospherics.  Measured  with  the  selective  amplifier  it  was  as 
follows: 

at  10  kc/sec  on  average  12  p.v  background  noise, 
atmospherics  ^  120  (peaks) 

at  3  kc/sec  on  average  10.8  pv  background  noise, 
atmospherics  s  20  pv^^.  (peaks) 

When  the  antenna  was  short-circuited  -  grid  of  the  first  cascode 
grounded  -  the  voltage  was  0.2  yve££  (at  Bu2).  The  antenna  was 

calibrated  by  the  method  described  in  the  1962  Annual  Report.  The 
following  relation  was  found: 

1  mv  induced  voltage  (Ue)  corresponds  to  B  =  1*75  •  10  (Weber/m2). 
Thus,  the  smallest  field  strength  measurable  with  the  built-in  instru¬ 
ment  is  B  .  =  8.75  *  10  ^ (rreber/mz  ) .  At  3  kc/sec  a  field  strength 

mm  i  rj 

of  B  =  2  •  lO--1-1  (Weber/m2)  can  be  measured  with  a  microvoltmeter  on  the 
basis  of  the  measured  background  noise.  This  is  possible  only  if  the 
noise  level  of  amtospherics  is  very  low  as  it  is  the  case  in  the  mine 
at  great  depths. 

Fig.  I  shows  the  field  strength  meter. 

Literature  Ll , 2 , 3 J  • 
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1,2  Magnetic  dipole  antennas  for  receivers 

Introduction: 


If  VLF  signals  are  to  be  received  in  solid  or  liquid  media,  the 
space  available  for  antennas  is  usually  very  limited.  In  the  present 
paper  the  problem  is  to  be  studied  what  antenna  type  is  suited  best  for 
recording  electromagnetic  V7aves  of  very  low  frequency  (f  %  15  kc/sec) 
in  solid  media,  i.e.  below  ground.  Another  limiting  condition  for 
the  antenna  size  is  the  fact  that  the  whole  receiver  shall  be  easily 
portable . 


1.  Comparison  of  antennas 


The  above-mentioned  conditions  show  that  the  electrical  and 
geometrical  length  of  an  antenna  in  any  case  will  be  very  small  as 
compared  to  the  wavelength  A.  of  the  signal.  The  smaller  the  antenna., 
however,  the  smaller  its  effective  height  and  its  radiation  resistance, 
i.e.  the  more  loosely  it  will  be  coupled  with  the  electromagnetic 
field  [l] .  The  problems  lie  in  the  fact  that  the  damping  Rv  by  the 
losses  of  the  total  antenna  circuit  is  to  be  kept  low  as  compared  to 
radiative  damping  caused  by  the  radiation  resistance  Rs,  this  means 
that  an  antenna  type  has  to  be  found  whose  efficiency  is  as  high  as 
possible . 

Theoretically,  an  electric  dipole  antenna  of  small  dimensions 

ought  to  be  sufficient  also  for  the  VLF  region  as  long  as  it  is 

tunable.  A  materializable  vertical  electric  dipole  having  an  effective 

height  of  2  m  for  f  =  10  kc/sec  has  a  radiation  resistance  of 
-6 

Re  ~  10”  Q.  The  induced  voltage  at  an  assumed  field  strength  of 
E  „„  =  50  mv/m  is  100  mv.  The  antenna  impedance  of  a  small  electric 
dipole  given  mainly  by  its  capacitive  resistance,  is  extremely  high 
and  resonant  tuning  is  possible  only  at  large  inductance.  The  useful 
amplifier  input  voltage  will  only  be  a  fraction  of  the  voltage  induced 
in  the  antenna.  By  means  of  a  corresponding  frame  antenna,  an 


optimum  adaptation  and  resonant  tuning  arc  easy  to  reach,  radiation 
resistance  and  the  effective  antenna  height,  however,  are  very  small. 
Electric  dipoles  for  VLF  signals  are  effective  only  in  the  form  of 
long  wire  antennas . 
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2.  Magnetic  dipole  ferrite  antenna 

The  introduction  of  a  magnetic  material  (ferrite)  makes  possible 
to  concentrate  the  flux  through  a  frame  in  a  red  of  much  smaller 


dimensions. 


Thus, 


eff.  height  h 


2mPMeff 

A 


(1,5) 


the  radiation  resistance 


Rs  - 


320AVu^ff 


(1,4) 


The  voltage  induced  in  the  antenna: 

Dind  ‘  cpS,m\|£j 


(1,5) 


When  timing  to  resonance  we  obtain 


Do  -  Dind'Q  -  2"fnpBpeffe. 


(1,6) 


Dimensioning  (Fig.  1,5) 

The  value  of  in  Eq.  (l,3)  depends  on  and  on  the  ratio  A/d 

in  Fig.  1,3  and  should  be  as  large  as  possible.  The  number  of  turns  n 
is  given  by  the  inductance  L  necessary  at  a  given  frequency  f.  Fig.  1,4 
shows  the  dependence  of  the  inductance  of  a  short  coil  (n  =  300 )  on  the 
one  hand  on  its  position  on  the  rod  and  on  the  other  hand  on  the 
rod  length.  Rod  length  A  =  32.5  cm,  curve  2}  A  =  65  cm,  curve  4$ 

A  =  97 • 5  cm,  curve  6. 

Detailed  studies  of  [4]  have  shown  that  it  is  favorable  to 
distribute  the  necessary  turns  over  a  winding  width  of  A/2  in  the 
rod  center.  Eq,  (l,6)  for  the  voltage  on  the  ferrite  antenna  is  then 
to  be  multiplied  with  the  reduction  factor  k,  since  the  winding 
concentrates  no  longer  at  the  rod  center,  having  the  permeability 
k  =  1  for  a  short  coil  in  the  rod  center 


k  =  0.7  for  a  long  coil  over  the  whole  rod. 

The  final  form  of  the  equation  for  the  voltage  on  a  tuned  ferrite 


antenna  then  reads: 


Uq  =  27ifnFBpeffQk 


(1,7) 
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b)  Design  of  an  optimum  ferrite  antenna 

The  core  material  1100  N22  of  the  Siemens  +  Halske  Co.  was 
chosen  for  the  reception  of  VLF  signals  in  the  frequency  region 
below  20  kc/sec.  Its  data  are  the  fallowing: 

>*tor  -  1300 
f  .  =1  kc/sec 

f  =20  kc/sec 

The  material  originally  intended  to  be  used  for  pot  cores  was 
pressed  into  tube  cores  for  antenna  experiments.  They  have  the 
following  dimensions: 

Length:  A  =  16  cm,  outer  diameter  d  =  19  mm 

inner  "  9  nun  . 

Optimum  rod  length: 

Fig.  1,4  shows  that  an  increase  in  rod  length  to  more  than  98  cm 
at  a  given  diameter  would  not  be  of  any  use.  The  optimum  number  of 
turns,  coil  shape  and  coil  position  on  the  core  were  determined  by 
experiment.  Fig.  1,5  shows  the  dependence  of  the  resonant  voltage  Ho 
on  the  number  of  individual  rods  (=  total  length)  -  curves  2,4*5 » 6  - 
and  on  the  number  of  coils  (each  having  n  =  500 )  arranged  around 
the  rod  center.  An  increase  in  rod  length  to  more  than  6  individual 
rods  and  an  increase  of  the  number  of  turns  to  more  than  6  coils 
(n  =  1800)  does  not  increase  the  resonant  voltage  Uo  considerably. 

The  quality  factor  Q  in  Eq.  (l,7)  depends  on  the  one  hand  on 
the  losses  in  the  oscillating  circuit  and  on  the  other  hand  on  the 
iron  losses  of  the  rod.  The  measurements  of  Q  as  a  function  of  the 
coil  position  on  the  ferrite  rod  reported  in  various  publications 
have  proved  to  be  incorrect  at  least  as  far  as  the  VLF  region  is 
concerned.  The  reasons  may  be  that  these  measurements  were  con¬ 
ducted  at  much  higher  frequencies,  f  >  500  kc/sec,  where  Q  is 
determined  almost  exclusively  by  the  iron  losses.  The  measuring 
device  for  studying  the  band  width  and  the  quality  factor  Q  is  shown 
in  Fig.  1,6.  The  RC  generator  G  adjustable  in  decades  is  loosely 
connected  with  the  resonant  circuit  LC  via  a  high-ohmic  resistance  R. 
The  resonance  voltage  is  applied  to  the  tube  voltmeter  RV  via  a 
cathode  follower  and  indicated  as  Uo.  Two  series  of  measurements 
were  made : 
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1.  Band  width  b  =  2Af  as  a  function  of  the  coil  position 
on  the  core,  number  of  turns  =  const,  f  =  const. 

As  L  decreases  towards  the  end  of  the  rod,  the  LC  ratio 
also  decreases  at  a  constant  number  of  turns.  Result: 
curve  bj  of  Fig.  1,7. 

2.  Band  v/ir;h  b  =  f(a).  At  a  constant  resonant  frequency  f, 
the  LC  ratio  was  also  kept  constant  by  increasing 

the  number  of  turns  towards  the  end  of  the  rod.  Result: 


curve  bg  of  Fig.  1,7* 

Measurement  of  b^:  At  a  given  frequency  of  10  kc/sec,  C  was 
tuned  to  resonance  and  the  generator  voltage  Ue  was  changed  so 
far  that  U0  =  100  mv.  Then  f  was  changed  by  f  until  UQ  was  70  mv. 
The  -./iicated  value  Af  is  b/2.  The  voltage  U  necessary  for  a 
constant  as  a  function  of  the  coil  position  (curve  b^)  is 
plotted  on  the  right  hand  side  of  Fig.  1,7.  The  result  shows  that 
the  band  width  and  Q  are  constant  over  2/3  of  the  rod  length, 
b  increases  only  towards  the  end  of  the  rod.  A  comparison  between 
b1  and  in  this  region  clearly  shows  the  effect  of  the  decreasing 
L/C  ratio.  Maximum  Q  at  approximately  2/3  of  the  rod  length  (as 
reported  by  a  number  of  authors  for  higher  frequencies)  could  not 
be  found.  An  increase  in  Q  at  the  expense  of  a  decrease  in  resonance 


voltage  was  not  observed  either.  A  maximum  number  of  turns 
distributed  over  %  a/3  left  and  right  of  the  rod  center  will  always 
favor  optimum  dimensioning  at  a  given  frequency  of  w  =  .  The 

coils  may  be  shifted  towards  the  rod  ends  only  in  case  the  desired 


band  width  is  to  be  larger.  Resistive  damping  of  the  antenna 


circuit,  however,  is  more  adequate. 


3.  Design  of  a  VLF  ferrite  antenna  for  the  frequency  range  of 
1-20  kc/sec 


In  view  of  rhe  test  measurements,  a  rod  (98  cm  long)  consisting 
of  6  tube  cores  (d  =  19  mm)  v/as  chosen.  The  front  planes  of  the 


individual  sections  were  ground  and  the  sections  were  pressed 
together  (without  adhesive)  by  an  axial  plexiglass  rod  with 
tightening  screws.  The  joints  between  the  cores  (y,x,x',y'  in 
Figs.  1,3  and  1,7)  do  not  affect  the  rod  quality. 

The  minimum  cut-off  frequency  is  to  be  1  kc/sec,  therefore  a 
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total  of  8  cross  windings  having  300  turns  each  were  attached;  at 
higher  frequencies  only  part  of  the  coils  are  used*  Measurement 
of  the  hand  width  of  the  antenna  circuit  with  the  same  device  as 
used  for  the  diagram  in  Fig.  1,5  yielded  the  following  values: 
f  =  3  kc/sec  ...  b  =  10  cps ,  fn  =  10  kc/sco  ...  b  =  75  cps. 

A  comparison  of  the  output  voltage  of  this  optimized  ferrite 

antenna  with  that  of  a  small  electric  dipole  at  equal  field 

strength  is  of  great  interest.  For  this  purpose,  Eq.  ( 1 , 7 )  is 

written:  „  „  „  , 

27t.n.F.E.)i  „f.Q.k 

U  = - 21± -  > 

0  A 


with  n  =  2400 

F  =  cross-sectional  area  of  the  rod  minus  cross-sectional  area 
of  the  bore 

F  =  2?2  .  10‘V  Q  =  f  /b  =  133 

E  =  50  mv/.  k  'v  0.8 

y.e££  ~  480  A  =  3  •  10^  m(=  10  kc/sec)  . 

Substituted  in  the  above  equation,  these  values  yield  an  output 
voltage  of  Uq  a*  260  mv  which  may  be  applied  to  the  high-ohmic 
.'.nput  of  an  amplifier.  This  voltage  is  much  higher  than  that  of  the 
short  antennas  described  above: 

Short  capacitive  dipole  -  100  mv,  but  on’y  a  fraction  can  be 

used  as  amplifier  input  voltage. 

Ironless  fr^me  antenna  1  m  diameter:  TJ,  .  ~  0.15  mv 

ind 

tuned:  U  ~  20  mv 


4.  Screening 


In  order  to  eliminate  electric  disturbances  having  a  capacitive 
effect  on  the  antenna  windings,  slotted  aluminum  or  copper  cylinders 
were  tested.  As  the  antenna  is  of  high  quality,  these  screens  had 
a  strongly  damping  effect  and  reduced  the  antenna  voltage  to  almost 
50$.  Screening  cylinders  causing  lower  damping  had  a  diameter  of 
no  less  than  40  cm.  Thus,  the  antenna  was  quite  unhandy  for  solving 
special  measuring  problems  in  mines.  A  slotted  screening  cylinder 


mo/lrt  A  ^  A  (IT  A  T  A^rA'W  A  -f*  A  V»  *J  MA  OA^A 

MkkUV  Vi  Lb  *  Vi  II  Vi  CL  *4  J.UUUXU  VUU  vvy^^i  UJ.XC  ”  C.WU 

turns  -  0.3  mn  CuL,  cylinder  diameter  15  cm  -yielded  much  better 
results .  The  wire  winding  was  cut  open  and  every  other  wire  was 
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soldered  to  a  collecting  bar  on  either  side.  The  principle  is 
shown  in  Pig.  1,8.  The  capacitor  thus  formed  yielded  a  frequenoy- 
dependent  screen,  damping  mainly  noise  of  high  frequency.  Damping 
at  a  much  smaller  cylinder  diameter  is  low. 

For  detailed  data  on  optimum  ferrite  antennas  see  our  VLP  re¬ 
port  [4,5]  . 

1.3  Transmitting  antonna  SAIX  at  St.  Certraudi 


The  transmitting  antenna  SAVIII  which  had  already  been 
described  in  the  last  annual  report  [2]  and  which  was  spanned 
out  in  a  vertical  plane  in  the  mine  of  St.  Gertraudi,  was 
replaced  by  an  improved  and  more  powerful  frame  antenna.  The 
antenna  area  (1600  mz),  howevar,  was  kept  constant,  but  the 
number  of  turns  was  increased  from  3  to  10.  This  means  an 
increase  in  magnetic  moment  to  more  than  the  threefold  as 
compared  to  the  previous  antenna. 

As  the  existing  galleries  and  shafts  had  to  be  used  for 
building  the  loops,  the  antenna  area  does  not  exactly  form 
a  plane.  However,  the  diagram  shows  that  the  deviations  of 
the  antenna  shape  from  a  vertical  plane  frame  are  not  very 
large.  For  distances  of  measurement  that  are  large  as  compared 
to  the  frame  dimensions,  this  antenna  acts  like  a  horizontal 
magnetic  dipole  whose  axis  and  the  NS  direction  form  an 
angle  of  approximately  40°. 

N 

X  X 
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In  order  to  reach  a  large  antenna  current  and  thus  also  a 
large  magnetic  moment  with  the  available  transmitter  whose 
output  power  is  1000  w,  the  resistance  of  the  total  antenna 
windings  had  to  be  not  higher  than  1  ohm.  In  this  case  an 
antenna  current  of  more  than  30  a  occurs  which  is  a  consider¬ 
able  increase  as  compared  to  the  former  antenna.  The  low 
antenna  resistance  was  reached  by  a  sufficiently  large  wire 
cross  section  (aluminium  v/ire  being  135  nun2  in  diameter). 

The  construction  of  the  antenna  loop  involved  several 
electric  and  mechanical  problems  that  had  to  be  solved  satis¬ 
factorily.  In  this  connection  the  problem  of  insulation  has 
to  be  mentioned,  as  it  was  decisive  for  the  actual  attainable 
antenna  quality  as  well  as  fer  the  attainable  resonant  current. 
Insulation  was  impaired  especially  by  the  fact  that  a  great 
number  of  insulators  were  required  in  the  narrow  and  windy 
galleries , which  because  of  the  high  atmospheric  humidity  are 
covered  with  a  thin  layer  of  condensed  water  mixed  with  dust. 

If  furthermore  it  is  taken  into  account  that  the  voltage  across 
a  turn  may  be  at  least  300  v  under  good  working  conditions, 
a  maximum  difference  of  approximately  3000  v  between  antenna 
and  earth  occurs.  The  occurring  surface  leakage  currents  de¬ 
teriorate  considerably  the  quality  of  the  antenna  circuit, 
thus  decreasing  the  antenna  current.  Pig,  II  shows  a  section 
of  the  antenna  construction;  it  can  be  seen  that  the  insulators 
are  large  so  that  sufficient  insulation  is  guaranteed  also 
under  the  poor  conditions  in  the  mine.  The  frame  is  tuned  t *> 
series  resonance  and  adapted  to  the  amplifier  output  in  a 
way  that  has  already  been  described  in  former  reports. 

Measurements  with  this  transmitting  antenna  were  conducted 
up  to  a  distance  of  7  km,  in  this  case,  however,  part  of  the 
distance  of  measurement  gees  through  air,  as  the  St.  Gertraudi 
mine  has  no  galleries  of  the  desired  dimensions.  These  test 
measurements  showed,  however,  that  it  is  favorable  to  place 
the  receiver  deep  below  the  earth's  surface  us  thus  the  at¬ 
mospheric  noise  level  is  attenuated  and  the  full  sensitivity 
of  the  receiver  can  be  utilized. 
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1.4  Direction-finder  antenna 


For  detailed  studies  on  the  direction  of  the  electromag¬ 
netic  field  in  the  neighborhood  of  the  transmitting  dipole, 
the  existing  receiving  antennas  [2]  have  not  been  suited. 
Preliminary  experiments  showed  that  the  deviations  in  field 
direction  caused  by  differing  rock  conductivity  (transition 
between  dolomite  and  lead  ore  of  better  conductivity)  are 
only  a  few  degrees.  Furthermore,  the  change  in  field  direction 
(AX)  which  is  due  to  a  rotation  of  the  transmitting  antenna 
(AS)  varies  according  to  the  conductivity  of  the  medium,  which 
also  made  a  more  suitable  receiving  antenna  desirable. 

The  new  direction-finder  antenna  was  constructed  such  that 
the  axis  of  the  receiving  antenna  rotates  freely  in  space  and 
can  furthermore  be  fixed  in  any  desired  position;  the  plane 
of  rotation  of  the  receiving  antenna  can  also  be  chosen  accord¬ 
ing  to  the  conditions  of  measurement.  The  necessity  of  choosing 
a  certain  plane  of  rotation  of  the  receiving  antenna  was 
evident  already  during  the  measurements  in  the  mines  of 
Bleiberg  [7]  where  the  directional  pattern  of  the  transmitting 
antenna  along  a  steep  ore  body  was  measured.  The  antenna  is 
set  up  in  such  a  way  that  the  base  plate  which  simultaneously 
serves  as  a  plane  of  reference,  is  fixed  horizontally  at  the 
beginning  of  a  measurement,  then  the  antenna  axis  is  adjusted 
in  north-south  direction  (Fig.  1.9).  This  position  serves 
as  the  starting  position  relative  to  which  a  field  direction 
or  a  change  in  field  direction  can  be  measured  with  an 
accuracy  of  +  1°. 

Fig.  1.9  shows  the  direction-finder  antenna  and  the 
arrangement  of  the  circular  scales  in  a  schematic  diagram. 

The  antenna  circuit  on  the  whole  agrees  with  that  of  former 
receiving  antennas.  The  ferrite  rod  is  approximately  50  cm 
long,  the  coils  being  arranged  at  its  center.  The  resulting 
resonance  elevation  of  the  receiving  voltage  can  be  fully 
utilized  by  adapting  the  high  resonance  resistance  of  the 
parallel  resonant  circuit  to  the  low-ohmic  input  of  the  meter 
by  means  of  a  suited  impedance  transformer.  The  circuit  of  this 


impedance  transformer  which  has  an  input  impedance  of  approxi¬ 
mately  10  megohms  is  shown  in  Pig.  1.10.  Such  a  high  input 
resistance  can  be  reached  by  the  simple  method  of  a  double 
collector  base  circuit  equipped  with  two  transistors  of  high 
current  amplification.  The  frequency  dependence  may  be  largely 
neglected,  since  the  circuit  3hows  a  constant  amplification 
factor  (v  =  0.95)  and  a  constant  input  impedance  (10  megohms) 
up  to  600  kc/sec.  The  current  (6  -  12  v,  1  mA)  ie  supplied  by 
the  indicator. 

The  impedance  transformer  and  the  resonance  tuner  together 
with  the  receiving  antenna  form  one  unit  so  that  the  connecting 
cables  can  be  made  as  short  as  possible.  For  measurements  in 
other  mines  where  there  is  no  powerful  stationary  transmitter 
available,  the  50  w  amplifier  described  previously  is  used.  As 
far  as  a  well-known  and  reproducible  angle  $  is  required  for 
such  measurements,  a  coil  whose  axis  can  be  adjusted  into  a  given 
direction  with  an  accuracy  of  approximately  +2°  is  used  as  a 
transmitting  antenna.  The  transmitting  antenna  is  aligned  in 
the  usual  manner  by  means  of  a  compass  and  mine  map.  The  angle 
is  based  on  1  =  0°  and  can  be  chosen  arbitrarily  in  order  to 
study  the  transmitting  antenna  or  the  dependence  of  the  field 
direction  (%)  on  the  alignment  of  the  transmitting  antenna 
(S  )  (see  2.5).  In  contrast  *-o  former  transmitting  antennas, 
this  new  design  makes  possible  to  rotate  the  antenna  axis 
also  in  vertical  direction.  This  represents  a  considerable 
improvement  as  numerous  measurements  frequently  have  to  be 
conducted  along  a  steep  ore  body,  i.e.  passing  from  a  higher 
to  a  lower  floor  jjfj  • 
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2.  Propagation  measurements 
Review 


In  our  third  annual  report,  measurements  in  several  mines 
and  different  types  of  rock  have  been  described.  The  studies 
conducted  in  a  pc+assium  and  salt  mine,  various  iron  ore  mines 
and  a  coal  mine  were  treated  in  detail.  The  results  obtained 
were  of  great  interest,  showing  that  the  methods  of  measure¬ 
ment  developed  in  the  Tyr  lean  mines  of  St.  Gertraudi,  Schwaz 
and  Lafatsch  are  applicable  also  to  other  types  of  rock, 
yielding  good  results.  Farthermore ,  the  obtained  results  of 
measurements  were  compared  with  the  existing  results ;  thus, 
facts  were  obtained  that  hold  for  a  general  theory  of  propa¬ 
gation  in  different  types  of  rock. 

The  most  important  points  that  result  from  the  measurements 
conducted  in  Northern  Germany  in  1963  shall  now  be  summari?;ed 
in  brief.  The  program  of  measurement  of  the  1964  report  is 
based  on  them. 

1.  The  measurements  were  conducted  in  rock  having  a  strong¬ 
ly  changing  conductivity  (Konrad  I  having  approximately 

10"^  mhos/m  and  Hansa  III  approximately  10  ®  mhos/m).  In 
these  cases  the  conditions  of  propagation  were  as  expected, 
in  the  second  case  a  very  strong  attenuation  occurred,  with 
transition  from  the  pure  near  field  to  the  far  field  being 
observed  at  20  m;  attenuation  in  the  second  case  was  very  small 
and  was  related  with  a  very  great  expansion  of  the  near  field 
which  extended  even  beyond  the  remotest  points  of  measurement. 
If  the  conductivities  are  so  small  and  the  frequencies  are 
below  10  kc/sec,  the  pure  near  field  of  the  transmitting 
antenna  may  be  several  kilometers. 

2.  Almost  all  measurements  were  made  in  very  inhomogeneous 
rock.  The  mineralization  on  the  whole  was  concentrated  to  a 
very  small  section  of  the  range  of  measurement  having  approxi¬ 
mately  the  shape  of  a  vertical  ore  layer  10  to  20  m  thick. 

As  the  conditions  of  measurement  as  well  as  the  geometrical 
conditions  we re  too  complicated  to  be  accurately  evaluable 
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by  theory,  a  great  number  of  measurements  had  to  be  made.  A 
comparison  with  the  results  obtained  under  homogeneous  condi¬ 
tions  enabled  us  to  determine  the  effect  of  inhomogeneities 
from  the  above  results.  For  this  method,  however,  very  compre¬ 
hensive  data  of  measurements  are  required  in  order  to  evalu¬ 
ate  unknown  conditions  of  propagation. 

3.  Measurements  of  the  received  field  strength  of  the 
GBR  transmitter  that  were  conducted  in  Tyrolean  mines  as  well 
as  in  mines  of  Northern  Germany  showed  different  attenuation 
in  horizontal  and  in  vertical  direction.  Attempts  of  explaining 
this  phenomenon  have  already  been  made  in  [2]  and  in  [6j ,  yet 
there  still  exist  doubts  as  to  the  correctnesL  of  these  ex¬ 
planations.  A  repetition  of  the  measurement  in  a  mine  was 
very  desirable,  where  field  strength  measurements  in  horizontal 
galleries  and  in  vertical  shafts  were  possible  beginning  at  the 
earth's  surface. 

2.1  Measurements  at  Bloiberg 

In  order  to  complete  the  above  program  of  measurements,  an 
excursion  was  made  into  the  lead  mines  of  Bleiberg  (Karnten). 
The  measurements  conducted  there  can  be  summarized  as  follows: 

1.  Propagation  measurements  through  homogeneous  or  weakly 
inhomogeneous  rock. 

2.  Measurement  of  the  directional  pattern  of  the  trans¬ 
mitting  antenna  and  determination  of  the  rock  con¬ 
ductivity.  These  measurements  also  made  possible  state¬ 
ments  on  the  near  field  and  on  the  far  field. 

3.  Measurement  of  the  decrease  in  field  strength  of  the 
GBR  transmitter  in  rock  in  horizontal  and  vertical 
direction . 

This  program  wa3  intended  to  serve  as  a  supplement  to  the  above 
studies  and  to  yield  explanations  of  unknown  phenomena  such 
as  the  difference  in  attenuation  of  GBR  signals  in  vertical 
and  horizontal  direction.  This  was  the  first  measurement  in 
a  large  and  strongly  mineralized  lead  mine. 
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ad  l)  The  propagation  measurements  wore  made  at  a  frequency 
of  3  kc/scc.  The  rock  lying  between  transmitter  and  receiver 
mainly  consisted  of  Yfetterstein  limestone,  at  individual  points, 
however,  mineralization  in  the  form  of  lead  sulfide  and 
sphalerite  was  found.  These  ore  inclusions  were  only  small  as 
compared  to  the  total  distance  of  measurement  (20  -  30  m 
mineralization  as  compared  to  400  -  800  m  of  measurement). 
Therefore  the  effect  of  this  mineralization  could  not  be  3een 
in  the  measured  curve.  Figur  2.1  shows  the  curve  of  the 
Bleiberg  measurements  (curve  1  and  Table  2.1). 


Table  2.1 


m 

J-0 

Ue  tvV3 

90 

85 

2400 

192 

20 

400 

250 

20 

190 

320 

15 

110 

385 

18 

63 

400 

8 

51 

440 

3 

36 

650 

3 

8 

860 

2 

_ _JL 

At  a  distance  of  approximately  400  m,  attenuation  proved  to 
increase  owing  to  a  change  in  rock.  For  reasons  of  comparison, 
a  curve  calculated  for  10~^  mhos/m  conductivity  was  plotted 
as  a  dashed  line.  This  shows  that  up  to  a  distance  of  400  m 
the  conductivity  was  smaller  than  10  ^  mhos/m,  whereas  at 
greater  distances  it  was  higher  than  the  assumed  value.  If  the 
dashed  curve  is  assumed  to  be  a  mean  value,  which  may  be  done 
for  approximation,  the  attenuation  factor  for  the  first  kilo¬ 
meter,  beginning  at  the  transmitting  antenna,  is  approximately 
66  db.  At  the  first  glance,  this  value  seems  to  be  rather  high, 
but  it  must  be  taken  into  account  that  this  value  only  holds 
for  the  immediate  near  field  of  the  transmitting  antenna.  The 
greater  the  distance  of  measurement  in  the  far  field  of  the 
transmitting  antenna,  the  lower  becomes  the  value  of  attenu¬ 
ation  per  kilometer.  This  shows  that  the  attenuation  depends 
on  the  distance  from  the  transmitting  antenna  and  it  is 
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only  in  the  far  field  that  attenuation  per  kilometer  becomes 
constant.  Yet,  the  above  value  is  quite  important,  as  will  be 
explained  ir.  detail  in  chapter  2.2  in  connection  with  the 
near  field  and  the  far  field.  The  curves  2,  2a  and  3  plotted 
in  Fig.  2.1  show  the  effect  of  the  different  attenuations  in 
the  far  field  and  in  the  near  field.  Curves  2  and  2a  give  the 
decrease  in  field  strength  of  the  GBR  transmitter  during  the 
horizontal  penetration  into  the  earth's  surface.  The  two 
curves  only  differ  by  the  fact  that  the  distance  from  the 
shaft  mouth  (m)  is  plotted  in  curve  2,  whereas  in  curve  2a  the 
depth  (m*),  i.e.  the  shortest  distance  between  the  site  of 
measurement  and  the  earth’s  surface  are  plotted.  Curve  3, 
however,  gives  the  decrease  in  field  strength  in  vertical 
direction. 


Table  2.2 

horizontal  measurement 


m 

in* 

UeTuVj 

560 

150 

64 

650 

175 

44 

750 

205 

40 

850 

240 

51 

920 

260 

28 

Table  2 .3 


vertical 

measurement 

m 

UeEuVj 

0 

400 

50 

70 

100 

50 

200 

15 

300 

15 

400 

12.5 

450 

11 

500 

7.6 

600 

5.6 

650 

_ LlO 

In  contract  to  the  . ^ o  ux  mv# uoux uo  xu  one  m xne s  ox 
Northern  Germany  [ 2 ]  ,  both  curves  show  approximately  equal 
attenuation.  By  means  of  extrapolation,  the  attenuation  of 
the  GBR  transmitter  (l6  kc/seo)  in  rock  having  a  conductivity 
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of  5 •10“^  mhos/m,  was  found  to  be  approximately  22  db.  In 
order  to  explain  this  discrepancy  between  the  values  of  atten¬ 
uation  obtained  from  propagation  measurements  with  the  trans¬ 
mitter  set  up  in  the  mine  and  those  obtained  from  measure¬ 
ments  of  the  GBR  transmitter,  we  may  state  the  fact  that  in  one 
case  the  attenuation  in  the  near  field  was  measured,  and  in  the 
other  case  in  the  far  field.  At  what  distance  the  far  field 
begins  and  in  how  far  the  distance  depends  on  the  conductivity 
of  the  medium  and  on  the  frequency  is  shown  on  the  one  hand  by 
measurements  of  the  radiation  pattern  of  the  transmitting  an¬ 
tenna  (chapter  2.4)  and  on  the  other  hand  by  theoretical  con¬ 
siderations  discussed  in  the  next  ohapter. 

Besides  the  studies  on  the  above  problem,  the  conducti¬ 
vity  of  rock  was  also  determined  in  the  mines  of  Bleiberg 
by  measuring  the  directional  pattern  of  the  transmitting  an¬ 
tenna  by  means  of  the  function  y(r,k).  Figur  2.2  shows  a 
graphic  representation  of  the  function  ^>(r,k)  for  several 
conductivity  values  and  for  a  frequency  of  5  kc/sec.  The  values 
of  measurement  in  Table  2.4  are  also  plotted  and  make  possible 
to  determine  the  electrical  conductivity.  Since,  however,  the 
different  points  of  measurement  were  made  in  rock  of  changing 
ore  content,  the  results  show  an  adequate  straying.  For  dead 
rock,  the  mean  value  of  conductivity  may  be  assumed  to  be 
10  ^  mhos/m,  whereas  for  more  strongly  mineralized  rock  this 
value  increases  somewhat,  being  approximately  5*10  ^mhos/m, 
as  was  shown  by  measurements  of  the  four-electrode  configuration. 
How  far  this  method  of  determining  y>(r,k)  makes  possible  to 
prove  spatially  confined  inhomogeneities  is  discussed  in 
detail  in  [7] . 

Table  2.4 

Frequency  of  measurement 
5  kc/sec 


m 

Wr,k) 

192 

1.55 

n  on 

V  #  V  V 

100 

1.80 

101 

1.97 

150 

1.93 

191 

1.11 
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2.2  Near  field  -  far  field 


For  a  better  understanding  of  the  experimental  results 
it  is  important  to  know  v/hether  a  measurement  is  conducted 
in  the  near  field  or  in  the  far  field  of  the  transmitting 
antenna.  The  conditions  of  propagation  differ  decisively 
for  these  two  cases.  V/hereas  the  near  field  is  characterized 
by  the  fact  that  the  component  of  the  electromagnetic 
field  prevails  as  compared  to  the  component ,  this  re¬ 
lation  of  components  i3  reverse  as  the  distance  of  measure¬ 
ment  increases,  i.e.  at  the  transition  from  a  pure  near  field 
to  the  far  field.  If  the  measurement  is  made  at  a  great 
distance  from  the  transmitter  or  more  precisely  at  a  distance 
from  the  transmitting  antenna  that  is  large  as  compared  to 
the  antenna  dimensions  and  the  wavelength*,  only  the 
$  component  of  the  magnetic  field  vector  is  left  over.  In 
propagation  measurements  in  the  high-frequency  region, 
this  case  alone  is  considered.  At  the  short  wavelength 
usually  studied,  the  effect  of  the  near  field  of  the  trans¬ 
mitting  antenr.^  may  be  largely  neglected,  as  it  is  too 
small  as  compared  to  the  total  range  of  measurement .  If  the 
propagation  conditions  of  VLF  waves  are  to  be  studied,  this 
otherwise  neglected  range  is  of  main  interest.  Considering 
that  the  wavelength  in  a  medium  of  average  conductivity 
(5  •10’’^  mhos/m,  i.e.  dolomite,  for  example),  is  approximately 
4  km  at  a  frequency  of  5  kc/sc-c,  the  near  field  extends 
over  a  range  of  at  least  10  km.  Then  the  transition  zone 
begins  which  at  distances  larger  than  that  goes  over  into 
the  actual  far  field.  An  increase  in  frequency  or  conducti¬ 
vity  causes  a  decrease  in  the  near  zone, which  must  not  be 
neglected,  not  even  in  this  case. 

A  principal  point  which  has  already  been  pointed  out 
(chapter  2.1)  is  the  fact  that  the  attenuation  of  an  electro¬ 
magnetic  wave  in  a  conducting  medium  is  much  larger  in  the 
near  field  than  in  the  far  field.  This  becomes  evident  from 
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the  expression  (2.1)  in  chapter  2. 5  for  the  magnetic  field 
strength  of  a  radiating  dipole  in  a  conducting  medium. 

At  small  distances  r,  the  first  terms  with  high  powers 
of  r  give  the  main  portion  of  the  field  strength.  At  very- 
large  distances,  however,  the  exponential  factor  is  de¬ 
cisive,  whereas  the  parenthetical  expression  assumes  a  con¬ 
stant  value  determined  by  the  conductivity. 

Furthermore,  the  equation  shows  that  at  0°,  which 
corresponds  to  the  component,  the  field  strength  decreases 
much  faster  towards  zero  than  at  =  90  .  In  order  to  illust¬ 
rate  this  behavior  the  magnetic  field  strength  |  Hi  was 
calculated  for  several  frequencies  beginning  at  1  kc/sec 
up  to  100  kc/sec.  Figs.  2.5  through  2,10  give  a  represent¬ 
ation  of  the  field  strength  trend  as  dependent  on  the  distance 
of  measurement  r.  The  curves  0°  and  J^=  90°  were  cal¬ 
culated  and  plotted  always  for  one  frequency.  Thus,  a 
comparison  of  the  H  and  components  is  possible  at  other¬ 
wise  constant  parameters.  It  can  be  seen  that  the  families 
of  curves  are  almost  consistent  even  at  larger  distances 
(5  -  10  km).  Especially  at  5  kc/sec  we  may  not  speak  of  a 
pure  far  field  at  distances  of  up  to  10  km.  Even  at  higher 
frequencies,  e.g.  100  kc/sec,  an  essential  quantitative 
difference  between  the  two  field  components  and 
becomaa  effective  only  at  a  distance  of  several  kilometers, 
which  corresponds  ho  a  distance  of  measurement  of  at  least 
10  wave  lengths.  At  the  first  glance,  this  result  may 
be  surprising,  but  it  can  easily  be  confirmed  by  experiment 
(  see  section  2.4). 

The  fact  that  the  examination  of  the  near  field  yields 
interesting  and  useful  results  shall  be  shown  by  means  of 
Fig*  2.10.  It  shows  how  the  attenuation  of  electromagnetic 
waves  depends  on  the  distance,  furthermore  it  shows  attenuation 
to  assume  a  comparatively  small,  constant  value  within  a  cer- 


rn«  n*A 
~  v*  *  *  O  w  • 


hiaj  u 


interpreted  by  the  fact  that 


a  certain  frequency  at  a  given  underground  path  of  trans^ 
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mission  yields  the  most  favorable  field  strength  at  the 
point  of  reception.  This  shows  that  not  every  frequency  may 
be  used,  but  that  the  applied  frequency  has  to  be  deter¬ 
mined  in  accordance  with  the  electrical  properties  of  the 
rock  and  also  with  the  distance  of  every  individual  measure¬ 
ment.  For  transmission  below  ground,  the  most  favorable 
frequency  of  transmittion  has  therefore  to  be  determined 
in  accordance  with  the  conditions  of  propagation. 

In  the  part  that  follows,  another  item  confirming  the 
difference  of  propagation  conditions  in  the  near  field  and 
in  the  far  field  shall  be  discussed.  In  chapter  2.1,  the 
rather  high  value  of  66  db  was  determined  for  "attenuation 
value/kra" ,  which  holds  for  the  immediate  near  zone  of  the 
transmitting  antenna.  In  contrast  to  this  value,  there  is 
that  of  22  db  for  the  attenuation  of  the  GBP.  transmitter 
measured  in  the  same  mine  under  equal  environmental  con¬ 
ditions,  despite  the  much  higher  frequency  (l6  kc/sec  as 
compared  to  3  kc/sec).  These  two  results  may  be  looked  upon 
as  being  reliable,  as  they  were  obtained  in  several  different 
mines  as  well  as  in  different  types  of  rock.  It  has  already 
been  mentioned  in  chapter  2.1  that  this  discrepancy  in  atten¬ 
uation  of  the  electromagnetic  wave  may  be  explained  by  the 
fact  that  one  measurement  was  made  in  the  near  field  and 
the  other  one  in  the  far  field  of  the  transmitting  antenna. 
Considering  expression  (2.1)  of  chapter  2.3  we  find  that  the 
terms  decisive  for  the  near  field  decrease  with  the  third 
and  fourth  power  of  r.  In  contrast  to  it,  attenuation  in 
the  far  field  is  determined  only  by  the  factor  e  /r. 


2.3  Measurement  of  the  direction  of  the  magnetic  field 
strength 

In  the  propagation  measurements  conducted  and  described 
so  far,  the  direction  of  field  strength  is  not  taken  into 
account,  only  the  amount  of  the  field  strength  being  measured. 


log  of  field  strength  IHI 


FIG.  2  A 


distance  (m) 


log  of  field  strength  IHl 


FIG.  2.5 


Parameter  <o 
Frequency  3  kcps 


log  of  field  strength  IHI 


FIG.  2. 7 


log  of  field  strength  IHl 


FIG.  2. 9 


FIG.  2.10 
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A  comprehensive  investigation  of  the  propagation  of  electro¬ 
magnetic  waves  in  the  rock,  however,  has  to  include  the  field 
direction  and  its  change  as  dependent  on  the  properties  of  the 
medium.  The  results  obtained  so  far  suggest  that  the  direc¬ 
tion  of  the  field  is  well  suited  to  give  an  explanation  of 
the  structure  of  the  medium  much  more  accurately  than  the 
absolute  value  of  the  field  strength.  Whereas  the  latter 
depends  on  the  sum  of  the  electric  properties  of  the  medium 
between  transmitter  and  receiver,  a  change  in  the  field 
direction  that  can  be  well  observed,  may  be  caused  by  an 
inhomogeneity  that  is  small  as  compared  to  the  entire  distance 
of  measurement.  The  conductivity  integrated  over  the  entire 
region  of  transmission,  and  thus  also  the  absolute  value  of 
the  field  strength,  do  not  change  noticeably  owing  to  such 
a  small  inhomogeneity . 

On  the  basis  of  the  equations 

IH|»^H*.+  hJl  =  m.h(r,$,k) 


and 

h(r,  j^k) 


-kr 
e _ 

r 


2k 

r 


2k2  ) 


+  sin 


'H 


1_ 

,4 


2k 

3 


(2.1) 


describing  the  magnetic  field  strength  vector  under  the 
assumption  of  an  unbounded  and  weakly  conducting  medium  Li], 
the  direction  (y)  of  the  field  vector  can  be  calculated 
from  the  two  components  and  Hjt .  In  general,  the  following 
relation  holds  for  the  angle  ^  (see  Pig.  2,il): 

■4 


.  Hjt 

t™  X  =  T* 


or  with  allowance  for  the  dependence  of  the  components 
and  on  the  alignment  of  the  transmitting  antenna,  i.e, 
on  the  angle  a?-,  the  relation 

.  JL 


tan  x 


h(>-0°) 


(2.2) 
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Substituting  the  function  f( r,k)  £7]  and  solving  the 
equation  (2.2)  for^,  we  obtain  the  desired  expression  for 
the  direction  of  the  field: 

,  tan  ,\ 

±Z  =  arctan  f(r,k)  *  (2-5) 

This  expression  shows  that  the  direction  of  the  field  at  a 
certain  point  having  the  distance  r  from  the  transmitting 
dipole,  depends  on  J' as  wo  13.  as  on  the  conductivity  of  the 
medium  in  between.  At  a  small  conductivity,  i.e,  if  the 
value  of  y>(r,k)  lies  between  O.5  and  2,  the  field  direction 
^  changes  almost  proportionally  to  the  direction  of  the  an¬ 
tenna  axis.  It  is  only  at  increasing  conductivity  that  is 
approximately  90°  in  a  wide  Grange.  This  case  corresponds 
to  the  conditions  of  propagation  characteristic  of  the  far 
field  of  a  transmitting  dipole.  For  the  parameter  y>(r,k) 
ranging  from  2  to  0.1,  the  expression  (2.3)  has  been  cal¬ 
culated  and  represented  in  F.lg,,  2.12.  The  Xield 
direction  measured  in  an  almost  homogeneous  rock  at  St. 
Gertraudi  is  well  consistent  with  the  curves  plotted  in 
Fig.  2.12.  The  direction-finder  antenna  described  in  chapter 
1.4,  which  allows  angular  measurements  to  bo  conducted  with 
an  accuracy  of  up  to  a  few  degrees,  was  used  as  a  receiver. 


FIG.  2.12 
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Table  2.5 
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2 . 4  Propagation  measurements  at  St.  Gertraudi 

The  newly  developed  direction-finder  antenna  was  used  for 
measurements  in  the  St.  Gertraudi  mine  in  connection  with  the 
problem  of  the  direction  determination  of  the  electromagnetic 
field  vector  H  in  a  homogeneous  or  inhomogeneous  rock, 
dealt  with  in  chapter  2.3*  It  is  the  aim  of  these  examinations 
to  make  statements  on  the  dependence  of  the  field  direction 
on  the  alignment  of  the  transmitting  antenna  and  at  the  same 
time  to  get  a  further  possibility  of  determining  the  near  field 
and  the  far  field. 

On  the  whole,  such  a  statement  is  always  based  on  the  pre¬ 
vailing  component  of  the  dipole  field.  This  co± responds  to 
the  determination  of  the  function  y(r,k)  introduced  in  i 7 J 
and  represented  in  Fig.  2.2.  If  the  determined  function  value 
is  smaller  than  0.1,  i.c.  if  the  component  is  10  times  as 
large  as  the  r  component,  we  may  consider  the  measurement  in 
the  far  field  to  be  of  good  approximation. 

This  method  has  been  expanded  by  allowing  for  the  size  of 
the  field  as  well  as  for  the  change  in  field  direction  at  the 
point  of  reception  as  dependent  on  the  direction  of  the  trans¬ 
mitting  antenna  axis.  The  derivation  in  chapter  2.3  shows  that 
the  field  direction  is  in  close  rolatic  .  to  the  electric  pro- 
uerties  of  the  medium  through  which  the  electromagnetic  *ave 
propagates . 
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Figures  2.13  and  2.14  show  the  results  of  measurement  obtained 
at  3  kc/sec  and  11  kc/sec  at  several  distances  from  the 
transmitting  antenna.  Throe  curves  are  plotted  in  each  polar 
diagram,  for  the  distances  450  m,  720  m  and  820  m.  The  scale 
along  the  periphery  gives  the  angle  of  measurement  In 
radial  direction,  the  receiving  voltage  was  plotted  in  a 
logarithmic  scale.  The  two  components  and  of  interest 
(see  Table  2.6)  were  obtained  from  the  four  points  of 
measurement  that  correspond  to  four  different  a?  values. 

Fig.  2.14  shows  the  predominance  of  the  component  over 
the  HJ7  component  even  at  a  measuring  frequency  of  11  kc/sec 
and  at  a  distance  of  800  n.  This  means  that  the  near  field 
reaches  far  beyond  this  range*.  Taking  into  account  that  the 
wavelength  in  the  St.  Gertraudi  dolomite  at  that  frequency 
is  at  least  1  to  2  km,  the  above  result  is  by  no  means 
surprising.  As  long  as  the  distance  of  measurement  is 
small  as  compared  to  the  wavelength'  in  rock,  we  cannot  speak 
of  a  far  zone. 


Table  2.6 

Voltage  of  reception  in  mv 
Frequency  3  kc/sec 


II 

450  m 

720  m 

820  a 

0° 

50.37 

7-54 

4-90 

K 

28.60 

6.80 

4,51 

60° 

21.50 

4-99 

3.60 

90° 

16.91 

_ wn 

5.05 

Frequency  11  kc/soc 

A?/r  = 

450  m 

720  m 

820  m 

0° 

14.55 

4.05 

2.25 

50° 

15.50 

5-55 

2.07 

60° 

10.90 

2.26 

1.61 

90° 

- 

..  1.45 _ 

■A-.??  . 

The  field  directions  which  belong  to  every  point  of  measurement 
and  which  fit  well  into  the  field  line  picture  of  a  dipole, 
have  been  plotted  in  the  two  figures.  The  table  2.5  in  chapter 
2.3  gives  a  comparison  of  the  direction  of  the  transmitting 
antenna  axis  with  the  corresponding  field  directions.  As  this 
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measurement  w as  conducted  in  the  near  field,  a  small  angle  of 
rotation  ^  of  the  field  corresponds  to  the  angle  of  rotation 
of  the  transmitting  antenna  axis  (see  also  Pig.  1.12  where 
the  values  measured  at  St.  Gertraudi  are  plotted  besides  the 
curves  derived  theoretically). 

In  the  foregoing  sections  it  ha3  been  pointed  out  that  a 
statement  on  the  homogeneity  of  rock  is  possible  by  measuring 
the  fi.„ld  direction.  Whereas  in  former  measurements  only  the 
quantity  of  the  field  strength  vector  determined  by  a  con¬ 
ductivity  that  was  averaged  over  the  entire  space,  had  been 
measured,  the  field  may  be  distorted  already  by  small  spatial 
innomogeneities .  The  measurements  of  St.  Gertraudi,  conducted 
in  an  almost  homogeneous  rock,  are  used  as  preliminary  studies 
for  a  detailed  investigation  of  the  above  problem.  The  next- 
step  will  be  similar  measurements  with  the  new  and  improved 
device  in  inhomogeneous  or  mineralized  rock.  In  this  connection, 
a  study  of  the  refraction  of  field  lines  along  the  rock  -  air 
interface  might  be  of  special  interest.  Such  a  study  might 
make  it  possible  that  the  propagation  path  of  VLP  waves  can 
be  followed  when  leaving  the  conducting  earth's  surface  and 
entering  the  atmosphere.  It  has  already  been  shown  that  a 
refraction  occurs  along  this  interface  and  that  the  field 
direction  there  is  not  consistent  with  the  purely  geometrical 
fi.ild  direction  that  is  expected.  Because  of  the  unsuited 
device,  detailed  measurements  so  far  have  not  been  possible. 

The  sites  of  measurement  at  St.  Gortrau-Ji  are  well  suited  for 
such  a  program  of  measurement,  as  it  is  possible  to  conduct 
measurements  in  field  direction  below  ground  as  well  as  along 
the  rock  -air  interface.  A  change  in  measuring  frequency  of  up 
to  100  kc/sec  would  also  ">ffer  the  possibility  of  studying 
the  behavior  at  this  conductivity  discontinuity  in  the  pure 
near  field  as  well  as  in  the  far  field.  It  is  assumed  that 
different  effects  will  occur  which  are  expected  to  be  of  great 
interest. 
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3*  Conductivity  measurements 

3.1  Direct  current  measurements 


In  continuation  of  the  studies  discussed  in  the  1963 
report,  conductivity  measurements  with  direct  current  w.re 
conducted  also  in  1964.  The  problem  from  what  distance  a  on¬ 
ward  the  effect  of  the  gallery  space  is  negligibly  small  when 
applying  the  Y7enner  method,  was  solved.  For  this  purpose, 
first  a  quantitative  estimation  regarding  the  required  distance 
between  the  electrodes  was  made  by  means  of  an  analytical  re¬ 
presentation  of  the  flow  pattern  (section  J.l.l),  secondly, 
the  effect  of  the  rock  inhomogeneity  on  the  result  of  measure¬ 
ment  was  studied  by  averaging  over  a  large  number  of  measure¬ 
ments  conducted  at  different  sites,  and  thus  the  result  of 
estimation  was  checked  by  experiment  (section  3.1*2).  A 
strictly  mathematical  solution  of  the  problem  was  found,  but 
it  has  not  been  numerically  evaluated  up  to  date  (section  3.1.3). 


3.1.1  In  Scientific  Report  Nr.  9  [8j  ,  an  expression  was  de¬ 
rived  giving  the  current  distribution  between  two  punctiform 
electrodes  in  a  homogeneous  full  space.  In  the  spherical  co¬ 
ordinates  r,  and  S'  it  reads; 


= 


j  = 


0  = 


■h  (- 


r  -  hcos 


ft- 


r  +  h  cosJ^'  \ 
g8™?  ==5  ) 


+  h^  -  2rh  cos  ^  \!r+h^  +  2rhcos-^^ 

0  (3.1.1) 


( 


h*sin 


s- 


..sin/J^ 


r»*) 


4n  ir2  +  h2  -  2rh  cos^*5  1r2  +  h2  +  2rhcoo  ^ 


From  it  we  calculate  the  portion  1^  of  the  total  current  I 
which  flows  through  a  circle  having  the  radius  R  and  the  center 
M,  whose  plane  is  perpendicular  to  the  line  connecting  the 
two  electrodes,  i.e.  ft-  ~  .  For  this  case  the  following  ex¬ 
pression  is  valid: 


3  ~  2 


i 


I*h 


2  2u  V(r*£  +  h2)^ 


For  1^  we  obtain: 


2u  R 


h  • 


"n  rdr  dy  =  l(l  -  — o  )  (3.1.2) 


^  2n>V  +  h‘  )3 


R2  +  h2 


I1.100 

Fig.  3.1*1  shows  the  quantity  b  =  - - —  which  is  the  portion 

of  the  current  (in  percents)  of  the  total  current  versus  the 
quantity  R/2h  which  is  the  ratio  between  the  radius  of  the 
circle  on  the  one  hand  and  the  distance  between  the  electrodes 
on  the  other.  For  b  =  10$  we  have  R/2h  =  0.25*  If  the  diameter 

of  the  gallery  (1.5  m)  is  substituted  for  R,  we  obtain  2h  =  6  m. 

At  that  distance  between  the  electrodes,  in  an  undisturbed 
space,  90$  of  the  current  flows  in  the  plane  given  by  n/2 

and  the  point  M  outside  a  circle  who.-se  radius  is  equal  to  the 

diameter  of  the  gallery.  For  a  justified  application  of  the 
full-space  expression  for  the  Wenner  method,  this  distance  of 
electrodes,  however,  is  not  sufficient,  as  the  probes  are  too 
close  to  the  electrodes  where  the  disturbing  effect  of  the 
gallery  is  larger  than  in  the  mentioned  plane.  In  section  3*1.2 
it  is  shown  that  the  distance  of  electrodes  for  the  Wenner  method 
has  to  be  at  least  20  m.  This  would  mean  that  in  an  undisturbed 
region  (according  to  R/2h  =  1.5/20  ^  0.075)  only  1$  of  the 
current  would  flow  through  the  circle  given  by  the  diameter  of 
the  gallery.  The  study  that  follows  below  shows  the  suitability 
of  the  Wenner  method  which  is  favorable,  because  the  relative 
error  of  voltage  measurement  becomes  small  owing  to  the  great 
distance  of  the  probes. 


3.1.2  The  voltage  profile  along  the  line  connecting  the  electrodes 
and  along  a  line  at  the  opposite  side  of  the  gallery  was  measured 
and  compared  with  the  theoretical  voltage  profile  for  the  cor¬ 
responding  conductivity.  The  measurement  was  made  such  that  one 
voltage  probe  was  fixed  at  M  during  the  period  of  measurement, 
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the  second  one  w as  fixed  to  the  measuring  points  in  between. 


H -  R  - 

E1  S1  -  M  S2  E2 


1 _  1 

5 - 

-  h  - 

- - 

-  h  - 

- > 

Using  the  denotations  of  the  above  Figure  (see  [2]  ,  1.10) 
the  following  relation  is  obtained  for  the  conductivity: 


6^  -T- 

0  2  7i U  ir 


R 


-  R£ 


(3.1.3) 


All  €  values  were  calculated  from  the  measured  voltage  values 
on  the  electrode  side,  and  from  them  the  arithmetic  mean  6  was 
calculated.  For  the  mean  value  5,  the  theoretical  curve  of  the 
voltage  trend  was  plotted  in  a  diagram. 

Results  of  this  measurement  in  the  Morgenschlag  (diameter 
1.5  ni)  are  represented  in  i'igs,  3.-»3>  3.1.5  and  3*1.6, 

where  A  is  the  calculated  voltage  tre.  1  for  6,  B  is  the  voltage 
trend  on  the  line  connecting  the  electrodes,  and  C  is  the 
voltage  trend  on  the  opposite  side. 

The  theoretical  curve  is  in  good  agreement  with  the  vol¬ 
tage  trend  on  the  connection  line.  The  correct  absolute  value 
of  conductivity  will  be  obtained  only  if  the  voltage  on  the 
opposite  side  of  the  gallery  is  also  equal  to  the  calculated 
voltage.  It  is  only  then  that  the  current  distribution  may  be 
assumed  to  correspond  to  that  of  the  homogeneous  full  space 
over  the  maximum  portion  of  the  space  considered.  For  2h  =  18 
or  24  m,  sufficient  agreement  is  reached,  whereas  for  2h  =  10  m, 
the  voltage  between  the  probes  at  the  electrode  side  is  still 
twice  the  value  of  the  voltage  on  the  opposite  side,  when 
applying  the  tenner  method  (^locations  of  the  probes  in  the  uia- 
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gram  are  marked  with  arrows). 

In  the  gallery  at  St.  Gertraudi,  whose  cross  section  is 
much  smaller,  sufficient  agreement  is  given  already  at 
2h  =  10  n  (Pig.  3.1.4). 

The  following  measurement  completes  the  obtained  picture. 

For  the  Wen ner  method,  the  ratio  Ug/Ug  of  the  voltage  between 
the  probes  on  the  electrode  side,  and  the  voltage  on  the  oppo¬ 
site  side  was  formed  for  different  distances  a  (electrode 
current  1=1  ma).  Pig.  3.1*7  gives  all  measured  ratios 
VUG  versus  the  distance  a  (Table  3-1.1) •  Furthermore,  the 
mean  value  calculated  from  the  measured  values  are  plotted. 

The  values  measured  for  small  a  show  considerable  deviations 
due  to  the  rock  inheraogeneity .  The  values  for  a  =  5  or  10  m 
were  measured  in  the  Wostschlag  whose  dimensions  correspond 
to  those  of  the  Morgenschlag.  Here,  the  inhomogeneity  apparent¬ 
ly  is  stronger  in  the  direction  perpendicular  to  the  gallery 
axis;  this  would  explain  the  difference  in  the  ratios  TJg/u ^ 
in  Table  3.1.1*  which  occurs  with  the  electrodes  once  being 
on  one  side  of  the  gallery  and  then  on  the  other  side.  This 
might  also  be  the  reason  for  the  comparatively  great  deviation 
from  unity. 

The  difference  between  the  large  and  the  small  gallery 
which  had  been  expected,  becomes  evident  by  the  measurement. 

The  deviation  of  the  measured  values,  however,  is  equal  for 
both  galleries. 

The  conductivity  measurements  with  direct  current  in  the 
mines  thus  yield  the  following  results:  All  sites  of  measure¬ 
ment  show  an  inhomogeneity  of  rock  that  exists  already  in  small 
regions.  Therefore  only  a  mean  value  of  conductivity  can  be 
measured.  For  determining  this  mean  value,  the  Wenner  method 
proved  to  be  well  suited  if  the  distance  a  is  sufficiently  large, 
under  the  present  conditions  >  7-8  m  (diameter  of  the  gallery 
1.5  m).  At  that  distance,  the  Wenner  method  is  of  sufficient 
accuracy  also  when  small  electrodes  are  used  which  are  easy 
to  fix. 
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Table  5.1  a 


Morgens chlag,  Westsohlag 


continued 


M 

a 

V 

V 

'G, 

s 

X 

V 

at 

1 

1 

\ 

2 

2 

10. 

■  5  7 

0.120 

0.091 

1.3 

0.108 

0.096 

1.1 

12 

8 

0.094 

0.081 

1.2 

0.095 

0.076 

1.25 

13. 

■  5  9 

0.080 

0.0?8 

1.0 

0.078 

0.070 

1.1 

15 

10 

0.117 

O.O90 

1.3 

0.110 

0.17 

0.9 

20 

10 

0.073 

0.055 

1.3 

0.077 

0.092 

0.8 

25 

10 

0.053 

0,042 

1.3 

0.052 

0.052 

1.0 

30 

10 

0.070 

O.O63 

1.1 

0.071 

O.O49 

1.4 

~  =  1.22 

35 

10 

0.133 

0.056 

2.4 

0.095 

0.105 

0.9 

UG 

40 

10 

0.185 

0.115 

1.6 

0.104 

O.I67 

0.6 

Table 

JLlJb 

Gertraudistollen 

K 

a 

U_ 

TL, 

U_ 

IT 

IL 

at 

E1 

Gi 

E 

E? 

E 

X 

"g 

d 

°G 

1.5 

1 

1.46 

0.18 

8.1 

1.59 

0.16 

9.9 

2.5 

1 

1.40 

0.15 

9.3 

1.75 

0.18 

9.7 

3.5 

1 

1.75 

0.31 

5-7 

2.14 

0.20 

10.6 

4.5 

-A. 

1.90 

0.26 

7.3 

1.67 

0.23 

7.3 

5.5 

1 

1.87 

0.31 

6.0 

1.87 

0.31 

6.0 

n 

uv 

6.5 

7.5 

1 

1 

1.76 

1.79 

0.14 

0.21 

12.6 

8.5 

1.36 

1.24 

0.22 

0.21 

6.2 

5.9 

7.1 

8.5 

1 

1.17 

0.31 

3.8 

1.48 

9.5 

1 

1.51 

0.39 

3.9 

1.30 

0.20 

6.5 

10.5 

1 

1.01 

0.244 

4.2 

1.06 

0.24 

4.4 

3 

2 

0.85 

0.40 

2.1 

0.72 

0.31 

2.3 

4 

2 

0.67 

0.42 

1.6 

0.96 

0.40 

2.4 

5 

2 

0.66 

0.37 

1.8 

0.82 

0.42 

2.0 

n 

6 

2 

0.91 

0.36 

2.5 

0.60 

0.35 

1.7 

1.96 

7 

2 

0.65 

0.22 

2.9 

0.60 

0.39 

1.5 

8 

2 

0.84 

0.53 

1.6 

0.66 

0.38 

1.7 

9 

2 

0.52 

0.36 

1.4 

0.62 

0.31 

2.0 

i 

| 

I 
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Table  3.1  b 
continued 


Gertraudistollen 


M 

at 

a 

\ 

\ 

uo 

\ 

UQ 

G2 

UB 

*0 

4-5 

3 

0.62 

0.58 

1.1 

0.58 

0.36 

1.6 

U. 

■C 

5.5 

3 

0.43 

0.28 

1.5 

0.46 

0.38 

1.2 

6.5 

3 

0.61 

0.36 

1.7 

0.40 

0.3  6 

1.1 

7.5 

3 

0.47 

0.27 

1.7 

0.40 

0.35 

1.1 

6 

4 

0.37 

0.35 

1.06 

0.30 

0.31 

0.97 

i 

3.1.3  Conductivity  measurements  by  the  potential  distribution  ' 

} 

of  steady  currents  ;j 

I 

I 

i 

The  current  distribution  of  spherical  electrodes  (homo-  ^ 

geneous  full  space)  or  semispherical  electrodes  (senispace)  i 

has  been  described  in[8j.  The  current  density  j  was  obtained  as  j 

a  gradient  of  the  potential  $  j 

•V  _ 

j  =  -grad  #  ,  , 

i 

which  can  be  found  by  elementary  considerations  and  which  must 
be  a  solution  to  the  equation  =  0  with  the  corresponding 
boundary  conditions.  In  [8] the  possibility  of  determining  the 
contact  resistance  between  the  electrodes  and  the  medium  has 
been  discussed  and  the  definition  and  determination  of  the 
penetration  depth  of  the  excitation  was  taken  into  account. 

The  conditions  that  must  be  fulfilled  for  applying  the  ex-  • 

pressions  to  the  homogeneous  semispace  have  also  been  dis-  , 

cussed.  The  expressions  considered  in  [8]|  thus  hold  for  a  plane 
interface  in  a  homogeneous  senispace.  Measurements  below  ground, 
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however, were  conducted  in  the  gallery  of  a  mine  and  it  has  been 
shown  by  experiment  that  a  measurable  excitation  occurs  on  the 
opposite  wall  at  least  at  large  electrode  distances.  How  far  the 
potential  lines  and  the  depth  of  penetration  are  affected  by 
the  curvature  of  the  gallery  wall  is  to  be  shown  in  further 
studies.  The  gallery  of  the  mine  may  approximately  be  looked 
tipon  as  a  o/linder  of  infinite  length  and  the  radius  =  a  .  The 
electrodes  are  to  have  a  distance  of  9  =  ?  from  the  cylinder 
axis  and  z  =  zq  from  the  z  =  0  plane  (Pig.  3.1.8). 


Without  the  cylinder  being  present,  their  potential  would  oe  given 
by: 

l/2  1  f ^2 

=  I/4tc [(l/(?i  +  (z  -  zoy)  -  i/(<^  +  (z  +  zq))  ](3.1.4) 


If  this  is  changed  so  that 


(5.1.5) 


where  $3  is  the  stray  field  conditioned  by  the  hollow  cylinder, 
then  &  is  a  solution  to  the  following  boundary  problem: 

A£>  =  0  (3.1.5a) 

(1)  -nxS 


& 

3?iQ=a 


.  _(ii>  w  +  aT) 

=n  ”  r  20  > 


3? 


?=a 


(3.1.5b) 
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thus , 


<j>=a 


(3.1.5c) 


(3.1.5c)  follows  from  the  fact  that  for  o  =  a  there  cannot 

■> 

exist  a  component  of  the  current  density  j  that  is  vertical 
to  the  cylinder  mantle, 

(3.1.5a)  expressed  in  cylinder  coordinates  reads 


a2  $  i  i  a!®  a2£> 

a?5”  9  as  Y  W  s?" 


o . 


A  particular  integral  of  this  relation  is  the  following  ex 
pression: 


Kc(mj>1)g(z) 


(3.1.6a) 


?1  “  W  +  9q  ~  2??0  cos  9 


m  ...  separation  parameter 


g(z)=  cos  m(z  -  z  )-  cos  m(z  +  z  ) 

K  (x)  ...  modified  Neumann 
o ' 

tion. 

It  will  be  seen  that  the  boundary  condition  %=  a  can  be  ful¬ 
filled  with  the  formulation  for  &S 


func- 


o 


o O  CO 

fr:(o)Kn(m?)ln(  ma)g(z)dm  cosn 

n=0  ~o 


(3.1.6b) 


(3.1.6c) 


I  (x)  ...  modified  Sessel  function  £  =  1,  £  =  2 

n  o  n 

r  ...  Neumann  number  n  =  1,2,... 

•n 

and  that  the  integrals  of  (3.1.6c)  are  convergent,  (3.1.6c) 
follows  from  (3.1.6b)  by  applying  the  addition  theorem 

CO 

K^m^)  =  »  (3.1. 6d) 

y>  a 
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For  <j>  =  a  the  following  expression  is  valid s 


oc 

=  -  J  £n(-l)n  J  l'n(m)K^(iaa)ln(ma)ing(z)dm*cosn  sp 


^  |V?1  +  (z  "  z0)2  +  (z  +  z0)2  £=a 

(except  for  the  constant  1/4tc 
which  is  omitted  here) 
oo  I 

L  £n  cosnjpJ  In(m^)Kn(mpo)g(z)dm  I  .  (3.1.7) 


For  mf  (m)  we  thus  obtain 
n  ' 


X  (m<r>  )l 1  (ma) 

mf  (n)  =  (-l)n  ?T7-A-T-7-T- 
nx  '  '  K'(ma)I  (ma) 

n'  '  nv  ' 


.  (3.1.7b) 


Hence  we  have 


£a « -l 


$■  )  £„  oosn  f 


Kn(m?0)r  (ma) 


^(ma, 


Kn(mf)&(z)dia 


(3.1.8a) 

This  means  that  the  total  field  ^5  is  given  by  the  expression 
<p=*  +  <J>S  = 


?i  -  h\ +  (z  +  z0) 

~  ..  .  rKn(ni>o^n(ma)  ..  , 


-Ji)  '“"S' 


TFSa)" - Kn(m?)g(a)am  . 


(3.1.8b) 


Choosing  the  integral  representation  for  the  first  part  of  (3.1.8b) 
like  that  of  (3.1.7) »  we  may  also  write 


§-*fr«n 


cosny  Kn(mj>n)lTi(m^)g(z)dm  - 


n 
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t  cosn 
n  J  . 


Kn(m^o)lA(na)Kn(in^) 


K'(ma 
n  ' 


g(z)  dm  V  = 


1  V,  „  n 

4n  L — cn  C0Sn^  J  K- 
n  {  n 


[ln(n>5-)K'(ma)  -  Ii(ma)Kn(n5)]g(Z)de 


(5.1.8c) 


>r  <jjb  we  thus  obtain  an  expansion  into  Fourier  series  in  cos 


$  -  i  IX  °°Sn  f  An 


(3.1.9) 


whose  coefficients  A  (^>5z)  are  given  by  the  integrals 

°r  K  (mo  )  r  _ 

A  (J»z)  =  j^T — “  |.In(n?)Kn(ma)  “  J-i!i(ma)Kn(ni^)Jg(z)dm  .  (5.1.9a) 

-l  n(aa)  "  j 


For  q=  a,  the  parenthetical  expression  under  the  integral  of 
(5.1.9a)  becomes  equal  to  che  Y/ronski  determinant  A (ln(ma) ,Kn(ma) ) } 
its  value  is  thus: 


A(ln(ma),Kn(ma))  =  -  l/ma  . 


(3.1.9b) 


This  is  the  problem  which  is  of  greatest  interest.  On  the  cylinder 
surface,  the  potential  is  given  by  the  following  expression: 


<f>( a»z,*p)  «  -  J. 


£n  cosnf  1  5 qra - rtm- 


(3.1.9c) 


The  convergence  of  the  integrals  in  (5.1.9c)  is  easy'to  prove. 


Writing 


ry%)«W 

An ( a ; z )  j  raK i (ma) 


dm  = 


i  n 


-  C'/o  dm  +  j  V-  dm  =  £  +  Jg,  (5.1.10) 
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the  estimation 


?  Kn(mfo)g(Z) 

K2I  -|J  mK' (ma)  dm 


0° 


< 


K  (mP  ) 
n  jo 

..  'mK1  (ma) 
M  nv  ' 


g(z) 


dm 


co 


K„H>> 


M  ****». 

n 


mK1 (ma) 


dm  (5.1.10a) 


is  obtained  for  J^;  if  M  is  sufficiently  large,  we  may  use  the 
asymptotic  expression  for  Kn  and  K^. 


K  (mP  )  J— 1—  e“  m?o 

n  \  ^  7tm 


K' (ma) 
n'  ' 


JX 

-)  NTcma 


-ma 


Thus, 


00 


J2  -S 


-m(yo  -  a) 


n: 


dm 


M 


N 


a 

So 


(5.1.10b) 


this  integral  being  oonvergent  fur  jj>o  4=  a*  As  the  integrand  in 
stays  finite,  i.e.  Cauchy's  principal  value  exists  for  n  =  0, 
convergence  is  guaranteed.  The  current  passing  through  a  cir¬ 
cular  plane  having  the  radius  R  in  the  plane  z  =  0,  is  a 
measure  for  the  depth  of  penetration.  This  current  is  given  by 

2k  R 


J(R>  -  dy  (3.1.U) 


0  a 
2k 


'/•St 


§  dS  df 

z=0 


A 

A 
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R  cP 


J(E)  -  2  '  Ii(»a)K0(mJ>)] 

•dmj)  dj)  .  (5.1.11a) 


m  sm  mz 

o 


a  0 


The  latter  expression  (3.1.11b)  is  valid  only  for  <  yQ.  For 
o  >  o  the  incident  field  'reads: 

■'  J  ®  lOO 


5^  =  cosn^>j  In(mj5o)Kn(mj))g(z)din 


CO 

/" 


yielding 


p  2k 

j(r)  -  -  j  f  ft?  a? 


R  2k 


a  0 
cP 


■  2^  f  K!(mf)'  [Ki(Ba)(?cIl(m?o)  -  aIl(ma))  - 


0 


-  I^(ma)(aK1(ma)  -  ^^(rn^)  ]  sin  mzodm 


2k 


oc 


0 


I0(,2i>0)KiCc“a)  ‘  Ko(lI%):Ci(ma) 


K' (ma) 
ov  ' 


^(mR)J  si: 


RK, (mR)  |  sin  mzodm  . 


Joh^o)  - 


(3.1.12) 


The  expressions  (3.1.12)  and  (3.1.9c)  can  now  be  evaluated 
numerically  which  is  planned  to  be  done  in  the  near  future. 


3*2  Measurements  with  samples 

3.2.1  The  frequency  dependence  of  £  and  e?  was  measured  by  means 
of  rock  samples.  The  measuring  arrangement  used  is  described  in [2], 
It  was  studied  whether  the  samples  have  a  high  surface  conducti¬ 
vity.  For  this  purpose,  a  guard  ring  capacitor  was  used.  Since, 
however,  vhe  results  with  and  without  guard  ring  are  consistent, 


I 
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the  rest  of  the  measurements  were  conducted  without  a  guard  ring 
as  balancing  the  guard  ring  circuit  is  very  time  consuming. 

The  results  of  the  measurement  were  evaluated  by  K.W. 

Wagner '  s  [9]  dispersion  theory.  In  this  theory  it  is  assumed 
that  the  frequency  dependence  of  £  and  6  is  determined  not  by 
a  single  relaxation  time  t,  but  by  a  relaxation  function  y(r). 
t  changes  continuously  from  0  to  «>,  whereas  ^(x)  is  determined 
by  the  distribution  function  k(x)  of  the  relaxation  times  * 

Y/agner  assumes  that  the  distribution  of  x  is  governed  by  a 
probability  function.  He  assumes  a  distribution  of  the  relaxation 
times  round  about  a  predominant  relaxation  time  xq,  the  density 
of  the  distribution  being  given  by  the  distribution  constant  b. 
Wagner's  distribution  function  has  the  following  form: 


k(x)dx 


»b  ~bz  z 
6 


(3.2.1) 


This  leads  to  a  relaxation  function ^>(x) 

*(t)  / 


-bV-  a  -  (^A0)e'Zd2 


(3.2.2) 


with  the  quantities  z  =  In  hit 

^  oo 

z  =  In  x/x 

o 

By  means  of  this  relaxation  function  we  obtain  the  following 

expression  for  the  real  part  and  for  the  imaginary  part  of  the 

x 

complex  dielectric  constant  £  : 


.  f.  k«b  -b2  z2  f  -bKn 

6  “  TT  e  °J 

o 


z  e  cosh(2b:>  z  -  l)n 


1 


£  - 


8  »k«b  ,2  2 
oo  -b  z 


V?t 


2  C  -b2  n2 
ojo 


cosh  n 
cosh  2b2  z  n 


-dnj 


cosh  a 


dn 


(5.2.3) 


with  n  =  z  +•  z  s  In  kt, 
o 
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These  expressions  contain  four  constants,  each  having  a 
physical  meaning:  b ,  tq  and  k.  is  the  dielectric  constant 

at  optical  frequencies,  b  gives,  as  has  already  been  said,  the 
density  >f  distribution  of  relaxation  times  about  the  pre¬ 
dominating  value  tq.  For  tq  the  following  expression  is  valid: 


To  2rcf 

m 


(5.2.4) 


with  f  being  the  frequency  at  which  £"  has  a  maximum,  k  is 
m 

a  measure  for  the  total  dispersion,  since 


6  =  £  (l  +  k)  or 

c  °°  ' 


k 


k 


(5.2.5) 


for  all  values  of  b  at  o>  =  0. 

A  good  reproduction  of  the  experimental  data  by  this  theory 
would  prove  that  the  measurement  was  not  affected  by  any  dis¬ 
turbing  fac+ors  such  as  for  example  those  that  may  occur  by 
electrode  polarization. 

For  examining  experimental  data  on  the  basis  of  \7agner's 
theory,  Y/.A.  Yager  presented  a  graphic  method  that  has  been 
applied  in  the  present  study  [lOj  . 


5.2.2  Graphical  method  by  \7.A.  Yager: 

A  new  variable  <f  defined  by  =  is  introduced,  with 

t'l  =  a"  at  f(a)  and  =  a"  at  f(x)» 

For  f  >  f  we  choose  f(a)  =  10  f(x),  then  we  have <  1. 
x  m  ' 

For  z  we  obtain 
x 

z  =  In  f  /f  . 
x  x'  m 

For  a  given  frequency  ratio,  cf  is  independent  of  and  k, 

but  is  a  function  of  b  and  z^.  Yager  calculated  cf  as  dependent 

on  b  for  various  values  of  z  .  From  two  experimental  <f  values 

we  may  thus  plot  two  curves  z  versus  b  yielding  a  point  of 

X1 

intersection',  b  and  z^  and  f^,  respectively,  are  determined  by 
this  point  of  intersection.  Furthermore ,  Yager  calculated  the 


3  -  16 


curves  p  "/&  k  and  (e’  -  £  )/£  k  for  values  of  b  <  2  and  for 
different  values  of  zq.  For  the  value  of  b,  which  has  already- 
been  determined,  the  curves  calculated  by  Yager  yield 
and  (£'  -  /f^k  for  different  zq  values.  For  equal  frequency 
values  f  we  then  plot  £" /&  k  versus  a"  and  (&'  -  £■  )/<£  k 
versus  <fgXp*  If  the  experimental  data  can  bo  represented  by 
Wagner's  theory,  cwo  parallel  lines  must  be  obtained  whose 
slope  is  given  by  £TOk.  The  point  of  intersection  of  the  ordinate 
with  the  straight  line  for  £,'  yields  the  value  of  <5^.  Thus,  all 
constants  of  Wagner's  theory  are  determined. 

3-2.3  Experimental  measurements: 

The  716-C  General  Radio  capacitance  bridge  was  used  for 
the  measurements  applying  a  substitution  method.  The  first 
balancing  was  made  with  the  grounded  plate  of  the  capacitor 
altogether  removed,  i.e.  only  the  intrinsic  capacitance  C 

8 

of  the  measuring  capacitor  being  connected.  According  to  the 
expression  given  in  the  General  Radio  716-C  Capacitance  Bridge 
Manual  we  obtain: 


Cp  . . .  constant  capacitance  in  the  second  bridge  branch.  For 
the  second  balancing,  the  capacitor  contained  the  sample: 

C  =  Cp  -  (c  +  C  ) . 

F  m  s 

According  to  expression  (5)  of  the  Manual  we  obtain 


C  =  AC  =  C'  -  C  =  C  , 
x  nr 


and  for  the  real  part  £,'  of  the  complex  dielectric  constant: 


Cq  was  calculated  from  the  geometrical  dimensions  of  the  sample. 


4*V»rt  1  Af«n 

x  wi.  vue  x vou 


1.1.  -  r.T  1  «  „  _ _ _ _  .*  _  _ 1  J  J  - 

uiic  wAjjiuooiuii  -l©  vailUt 


,<*>)•.  I  b 
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D 

x 


D') 


C' 


D, 


(3.2.7) 


as  D'  was  zero  in  all  measurements. 

The  complex  dielectric  constant  has  the  form 


the  loss  factor 


6*  =  a'  -  ife“, 

tan  <f  =  D  =  ~  . 
x 


If  (3.2.6)  is  substituted, 


n 


(3-2.8) 


is  obtained  for  the  imaginary  part  of  the  complex  dielectric 
constant.  As  the  curves  calculated  by  Yager  comprise  only  a 
certain  range  of  frequencies,  which  depends  on  f  ,  not  all  curves 
could  be  evaluated.  Especially  the  humid  samples  that  are  most 
similar  to  the  natural  state  of  the  rock  are  unsuited  for  an 
evaluation;  their  D  values  are  very  high  at  low  frequencies, 
therefore  they  cannot  be  measured  with  the  available  bridge. 

Very  dry  samples  for  which  f^  is  much  smaller  than  1  cps  are 
also  unsuited  for  evaluation,  as  the  measuring  bridge  cannot 
be  used  fox’  f  <  30  cps.  For  medium  values  of  humidity,  however, 
all  four  constants  of  'Vagner's  theory  could  be  determined. 

In  the  present  report  the  result  of  measurement  obtained 
with  a  dolomite  sample  that  was  studied  at  seven  different  values 
of  humidity,  shall  be  discussed.  Humidity  was  changed  in  the 
same  way  as  described  in  [jLlJ.  In  [ll]  it  was  stated  that 
especially  the  humid  samples  showed  certain  difficulties  in 
reproducing  the  measured  values  as  well  as  a  considerable  change 
of  the  measured  value  during  the  measurement.  This  could  be 
avoided  by  storing  the  samples  for  several  days  in  a  room  of 
constant  absolute  atmospheric  moisture.  In  order  to  increase  the 
moisture  content  of  the  samples,  the  absolute  atmospheric 
moisture  was  increased  by  raising  the  temperature  ,  the  relative 
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atmospheric  moisture  being  kept  at  100$.  Attention  must  be 
paid  to  the  fact  that  the  samples  had  to  be  stored  under  constant 
conditions  for  a  sufficient  period  of  time  in  order  to  keep 
them  uniform  as  to  their  moisture  content. 

The  sample  discussed  here  is  a  dolomite  sample  of  low 
porosity.  It  was  first  dried  (measurement  l)  and  then  stored 
in  a  humid  atmosphere  (measurements  2-6).  Then  the  sample 
was  kept  under  water  for  one  week  (measurement  7)*  Keeping 
the  sample  under  water  for  a  longer  period  did  not  increase 
its  water  content.  The  amount  of  contained  water  was  determined 
by  weighing,  it  is  expressed  by  the  value  p  =  ■ 

0  =  =  weight  of  the  dry  sample 

AG  =  G„  -  G. 

N 


3.2.4  Results  of  measurement 

The  results  of  the  seven  measurements  of  the  sample 
(dolomite)  at  different  values  of  water  content  are  shown  in 
Figs.  3*2.1,  3*2.2  and  3*2. 3j  respectively,  and  C$  being 

plotted  versus  the  frequency  f .  The  values  for  are  obtained 
from  those  for  <c"  calculated  from  the  complex  dielectric 
constant. 

£*  =  f0(£'  -  06")  -  e  (e.'  -  jjf-)  . 

O 

This  yields: 

£"  *<£0*a>  (3*2.9) 

The  S'  values  were  plotted  to  make  possible  a  comparison  with  the 
measured  values  of  [2l  .  At  all  humidity  values,  S'  was  found 
to  increase  with  frequency,  but  the  increase  becomes  smaller  as 
the  water  content  increases.  This  fact  might  be  an  explanation 
of  the  observation  that  the  conductivity  measurements  made  in 
solid  rock  (section  3*3)  and  the  measurements  reported  in  ^12^ 
and  do  not  remarkably  depend  on  the  frequency,  since  it  is 

not  known  whether  the  water  content  of  the  sample  after  the 
treatment  required  for  its  production  can  be  restored  to  its 
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original  value.  At  a  further  increase  in  humidity,  an  even 

smaller  increase  of  ^  with  frequency  is  possible. 

All  four  constants  of  Wagner's  theory  could  be  determined  only 

for  the  measurements  Nr.  3»  4>  5  and  6.  For  measurement  Nr.  1 

(dry  sample),  only  the  value  of  b  could  be  estimated,  whereas 

for  Nr.  2  an  estimation  was  possible  also  for  f  .  The  estimation 

r  m 

given  for  measurement  Nr.  7  (after  saturation  in  water)  is  very 
unexact,  since  the  diagram  b  versus  z  gives  a  grazing  inter¬ 
section. 

In  3«2.2  it  has  already  been  mentioned  that  the  straight¬ 
ness  of  the  curves  £"/<£  k  versus  and  (/='  -  £  )/s  k  versus 

'  oo  exp  °  00' '  00 

^exp  a  measure  the  accuracy  with  which  Wagner's  equations 
explain  the  experimental  data.  For  measurements  Nr.  3,  4>  5  and 
6  the  respective  diagrams  are  shown  in  Fig.  3*2.4.  The  points 
of  measurement  al  .  lie  on  a  straight  line  with  sufficient 
accuracy  to  make  possible  the  determination  of  <5^  and  k.  The 
results  are  summarized  in  Table  3*2.1.  The  first  column  gives 
the  number  N  of  the  measurement,  the  second  one  the  weight  G 
of  the  sample,  the  third  one  the  increase  in  weight  P  in 
percents  as  compared  to  the  weight  G^  of  the  dry  sample. 


Table  3.2.1 


N 

G  [g] 

m 

b 

fffl(cps) 

4o 

k 

1 

34*9518 

0 

a*  0.10 

2 

34*9523 

0.00143 

0.12 

0.017 

3 

34*9531 

0.00372 

0.19 

0.36 

9.2 

27.8 

4 

34*9560 

0.0120 

0.24 

11 

11 

24*4 

5 

34.9563 

0.0129 

0.26 

41 

16 

14.7 

6 

34.9569 

0.0146 

0.29 

81 

17 

11.8 

7 

35.0157 

0.183 

>  1.0 

%  800cps 

The  table  shows  that  the  four  constants  of  Wagner's  theory  depend 
on  the  water  content  of  the  rook.  The  •value  b  e-ivin./r  the  densit-v 
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of  the  distribution  of  relaxation  times  about  the  predominating 
value  increases  with  the  moisture  content.  This  corresponds 
to  a  decrease  of  the  frequency  range  of  anomalous  dispersion. 

The  value  f^,  i.e.  the  frequency  at  which  has  a  maximum, 
also  increases  with  the  moisture  content.  Only  the  value  k 
which  is  a  measure  of  the  total  dispersion  between  f  =  0  and 
f  =  oo,  decreases.  The  increase  in  £ '■  with  p  can  be  explained 
by  the  mixed  body  theories,  since  the  value  of  the  dielectric 
constant  increases  with  the  portion  by  volume  of  the  constituent 
that  has  the  highest  dielectric  constant,  in  this  case  it  is 
water  having  3=  81. 

Pig.  3*2.5  shows  the  frequency  dependence  of  the  experimental 

values  e*  and  as  well  as  the  values  and  <g Ji.  of 

"exp  ‘"exp  W  w 

measurement  6,  calculated  by  Wagner's  theory,  since  the  curves 
oalculated  by  Yager  comprise  best  the  used  frequency  range 
for  this  measurement,  whereas  for  the  other  measurements  they 
comprise  only  smaller  parts  of  the  frequency  range.  In  order 
to  illustrate  the  effect  of  the  constant  b  on  the  result,  the 
data  for  the  correct  value  of  b  =  0.29  (solid  line)  and  for 
the  arbitrarily  increased  value  b  =  0.31  (dashed  line)  were 
plotted.  The  measured  points  of  the  two  quantities  £'  and 
are  well  reproduced  by  the  curves  calculated  with  b  =  0.29. 

Too  high  a  value  of  b  yields  a  steeper  slope  of  the  curves 
and  thus  also  a  frequency  dispersion  of  6 which  decreases 
with  the  moisture  content,  since  is  calculated  by  multiplying 

with  u)£  . 
o 

This  examination  was  also  made  with  other  samples  of  rock, 

but  only  for  2  values  of  humidity.  A  measurement  conducted  with 

a  second  dolomite  sample  yielded  the  same  results  as  above.  Also 

for  new  red  sandstone  used  for  studying  the  reproducibility  [11], 

the  constants  b,  f  ,  £  and  k  assume  similar  values  as  for 

m  00 

dolomite,  as  well  as  the  samples  described  in  [2]  from  the  mine 
at  Willroth,  which  consist  of  ore  embedded  in  dead  rock.  All 
these  different  sorts  of  rock  have  approximately  the  same  po~ 
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rosity.  The  results  obtained  for  a  sample  from  the  ore  mine 
of  Salagitter  (oolitic  ore)  cannot  be  explained  by  Wagner's 
theory.  In  Fig.  3*2.6  they  are  shown  for  t wo  different  values 
of  water  content. 

In  contrast  to  the  other  samples,  the  curves  fcr  show 
a  minimum  that  cannot  be  explained  by  Wagner's  theory.  The 
dispersion  cannot  be  represented  by  a  relaxation  function  for 
which  a  distribution  of  relaxation  times  about  a  prevailing 
value  of  is  assumed.  The  shape  of  the  curves,  however,  would 
rather  suggest  two  different  prevailing  relaxation  times.  An 
evaluation  in  this  sense  has  not  been  attempted  so  far. 

'/  On  the  whole  it  may  be  said  that  the  frequency  dispersion 
of  conductivity  and  dielectric  constant  of  mineral's  observed 
in  rock  samples  can  be  well  described  by  Y/agnor '  j  dispersion 
theory.  As  this  theory  only  contains  four  constants  which  all 
have  a  physical  meaning,  and  as  it  is  based  on  the  Curie  law 
saying  that  the  reversible  absorption  current  is  proportional 
to  the  applied  voltage  and  to  the  geometric  capacitance,  this 
agreement  may  be  taken  as  a  proof  of  the  fact  that  the  sample 
measurement  yields  correct  values  for  and  and  thus  also 
for  6/,  and  that  the  results  of  measurement  are  not  affected 
by  other  effects. 


3.3  Measurements  with  alternating  current  in  solid  rock 


3.3*1  The  716-C  General  Radio  Capacitance  Bridge  was 
also  used  for  measuring  the  propagation  resistance  between  twe 
electrodes  fixed  to  the  rock,  namely  steel  pins  described  in  [2^ 
which  have  also  been  used  for  the  d-c  measurements  described  in 
section  3*1»  They  have  a  diameter  of  2a  =  0.7  cm  and  a  length 
of  h  =  5  cm.  The  propagation  resistance  of  such  an  electrode 
can  be  calculated  easily  (see  e.g.  P8p-  We  obtain 


R  = 


2ixh  6 


In 


2h 

a 


(33*1) 


From  a  certain  electrode  spacing  onward,  tr.o  total  resistance  of 


I" 


ste  ou 


100  f  l kcpsJ 
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the  arrangement  is  practically  independent  of  that  distance. 

This  is  the  case  when  the  electrical  flow  pattern  near  one 
electrode  does  not  change  when  the  electrode  spacing  is  varied. 
The  total  resistance  then  is  twice  that  of  one  electrode 

R  =  ^hlnf-  <5.3.2) 

The  above  expression  is  the  same  as  that  used  in  [l3] .  The 
constant  factor  which  is  contained  therein  and  which  causes  a 
decrease  in  resistance,  is  due  to  the  fact  that  for  the  used 
dipole  arrangement  the  distance  between  the  two  dipole  halves 
still  has  a  certain  effect  on  the  resistance. 

A  Schering  capacitance  bridge,  e.g.  the  716-C  General 
Radio  Bridge  can  be  used  for  determining  the  quantities  C  and  R 
of  an  equivalent  circuit  diagram  of  a  dissipative  capacitor.  As 
the  capacitance  between  two  electrodes  of  the  above  form  is 
very  low  for  rock  with  a  small  dielectric  constant,  parasitic 
capacitance  of  the  feeder  cables  and  bridge  screening  are 
usually  higher  than  the  capacitance  between  the  electrodes. 
Their  elimination  by  a  substitution  method,  however,  is 
difficult  when  measuring  in  solid  rock,  since  the  measuring 
bridge  during  the  measurement  is  located  opposite  a  wall 
having  not  a  uniform  potential,  but  a  potential  drop  along  it. 
This  makes  impossible  to  get  an  accurate  measurement  of  the 
capacitance  between  the  electrodes.  The  quantity  of  the  re¬ 
sistance  R  of  the  equivalent  circuit  diagram,  however,  can  be 
determined  accurately,  because  it  is  not  affected  by  the 
geometrical  arrangement  of  the  feeder  cables  in  contrast  to 
the  total  capacitance,  i.e.  the  sum  of  the  parasitic  capacitance 
and  the  capacitance  between  the  electrodes.  This  can  also  be 
shown  by  means  of  measurements,  e.g.  by  conducting  two  measure¬ 
ments  with  the  same  geometrical  arrangement,  where  one  electrode 
first  is  located  on  the  potential  of  the  bridge  screen,  then  it 
is  used  as  a  hot  electrode.  V/hereas  the  value  of  capacitance 
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differs  considerably  in  the  two  measurements,  the  value  of  the 
resistance  R  remains  constant. 

3.3*2  Results: 

In  contrast  to  the  sample  measurements  (section  3*2), 
the  results  of  conductivity  measurements  in  solid  rock  do  not 
yield  a  remarkable  frequency  dispersion.  They  are  shown  in 
Pig.  3*3.1 .  The  values  are  calculated  from 

<5-5-3> 

Because  of  tho  limited  resistance  range  of  the  capacitance 
bridge  at  low  frequencies,  d  could  only  be  determined  between 
5  and  100  kc/sec.  Pig.  3.3*1  shows  five  measurements  made 
at  electrode  spacings  of  d  =  0,5,  1,  2,  3  and  4  m.  The 
dependence  of  conductivity  on  frequency  is  very  small  and  the 
value  of  d  is  of  the  same  order  of  magnitude  as  the  d-c  value 
determined  by  a  four-electrode  arrangement.  The  small  fre¬ 
quency  dispersion  and  the  resulting  absolute  value  are  in 
good  agreement  with  the  results  of  conductivity  measurements 
in  rock  described  in  [l3j, which  were  conducted  in  drill  holes. 

In  the  above  report,  H  was  not  directly  measured  with  a  bridge 
because  of  the  long  connection  wires,  but  the  observed  impedance 
was  substituted  by  an  equivalent  circuit.  From  it,  R  was  de¬ 
termined.  This  complicated  method  is  superfluous  in  the  method 
described  in  the  present  report,  since  it  does  not  require  long 
cables  causing  a  considerable  loss  factor?  for  the  present 
method  bare  wires  wore  used,  but  the  principle  of  measurement 
was  equal  to  that  of  the  above  report. 

The  results,  however,  do  not  agree  sufficiently  with  the 
d  values  measured  in  samples,  neither  with  respect  to  the  ab¬ 
solute  amount  of  d,  nor  to  the  frequency  dependence.  The  reason 
for  the  differences  in  the  results  of  the  two  methods  of 
measurement  has  not  been  explained  satisfactorily,.  Above  all 
it  seems  necessary  to  conduct  measurements  with  humid  samples 
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and  to  extend  the  measurements  in  solid  rock  up  to  the  fre¬ 
quency  region  of  30  cps  -  200  kc/sec,  by  improving  the 
measuring  device,  in  order  to  conduct  an  accurate  comparison 
of  the  results.  This  shall  be  done  in  continuation  of  the 
present  studies.  Furthermore,  an  attempt  shall  be  made  to 
apply  the  evaluation  of  measurements  by  Wagner's  theory  also  to 
rock  samples  having  their  natural  water  content. 


4-1 


4  •  Theoi'etical  part 
4,1  Model  tests 


The  great  interest  in  the  propagation  and  excitation  of  electro¬ 
magnetic  waves  of  low  frequencies  as  well  as  the  technological  and 
financial  difficulties  involved  by  VLF-wave  studies  favored  the 
use  of  model  tests.  Such  tests  have  already  been  conducted,  e.g.  infv-i] 
Some  general  considerations  on  the  results  of  model  tests  in  conducting 
material  of  this  special  type,  which  are  intended  to  describe  the 
propagation  below  ground,  however,  show  that  they  are  not  of  very 
great  importance* 

Discussions  on  such  model  tests  have  already  been  published  injl'Sj 
and  shall  be  repeated  in  a  more  critical  form  in  the  present  paper. 

The  wave  equation  for  the  dimensionless  direction  vector  yields 
the  relations  between  the  scale  factor  oc  ,  the  electric  parameters  £ , 
yu,  and  G  of  the  simulating  model  and  the  model  medium,  as  well  as 
the  frequency 


<*=  r/rm  =  (tjt)1'2  (4,1a) 

6js  =  (4,1b) 

The  index  m  indicates  the  model  size,. 

On  the  whole,  only  one  q\;antity,  namely  G  of  the  equations 

in 

(4,1a)  and  (4,1b)  can  be  chosen.  It  determines  above  all  the  maximum 
attainable  model  factor ct , being  itself  determined  by  the  maximum 
attainable  model  conductivity  (in  such  a  way  that  the  medium  is 
in  a  state  in  which  tests  below  ground  are  still  possible). 

Assume  that  the  most  favorable  model  material  representing  a  hcmogen- 
ecus  full  space  hso  been  found  and  a  model  factor  oO  can  be  reached  which 
makes  possible  measurements  over  a  distance  of  up  to  100  km  in 
laboratory  tests,  there  are  still  several  difficulties.  In  this 
case,  antennas  whose  dimensions  are  much  larger  than  the  used  wave 
length  would  be  required  for  studying  the  process  of  propagation. 

The  model  test  would  then  yield  results  on  a  process  which  could 
never  be  materialized,  as  no  antenna  can  be  built  whose  length  is 
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a  multiple  of  the  wave  length.  If  the  conditions  were  so  favorable 
that  the  antenna  dimensions  are  only  fractions  of  the  wave  length 
also  in  the  model  test,  the  emission  of  such  antennas  would  be 
proportional  to  their  length  and  the  excitation  would  decrease  so 
rapidly  that  also  in  the  model  te3t  only  the  near  field  could  be 
studied.  In  situ  such  studies  do  not  involve  any  additional  technolo¬ 
gical  and  financial  expenses.  All  what  has  been  said  above  regards 
the  homogeneous  full  space. 

A  simulation  of  the  homogeneous  full  space,  i.e.,  investigation 
of  VLF  propagation  below  ground  and  at  a  great  distance  from  the 
earth's  surface,  in  the  optimum  case  would  only  represent  a  measure 
of  estimating  the  underground  homogeneity.  Although  this  is  of 
great  interest,  it  does  not  justify  the  great  expenses. 

It  would  be  much  more  important,  however,  to  simulate  the 
system  air  -  earth,  since  the  conditions  at  the  earth  -  air  inter¬ 
face  which  are  of  great  practical  importance  would  thus  be  made 
accessible  to  detailed  studies.  Here  too,  the  model  factor  <36  has  to 
be  chosen  in  such  a  way  that  the  conductivity  of  the  model  ground 
does  not  exceed  the  conductivity  of  water,  acidified  water  or  any 
other  liquid  exhibiting  favorable  properties.  Otherwise,  underground 
measurements  would  be  impossible.  The  possibility  of  simulating  the 
ground  by  metals  and  measurements  in  boreholes  has  already  been 
mentioned  in  [4  5 ~j  >  at  the  same  time,  hov/ever,  we  have  pointed  out 
that  the  measurement  can  no  longer  be  adapted  to  reality  owing 
to  secondary  effects,  e.g.  skin  effect.  The  use  of  a  comparatively 
small  maximum  conductivity  (of  liquids)  requires  the  choice  of  a 
low  model  factor  ct  which  means  that  only  the  near  field  can  be  studied 
if  the  dimensions  of  the  model  are  to  be  kept  in  tolerate  boundaries. 
The  latter,  however,  is  necessary  as  otherwise  the  model  test  would 
in  no  way  be  superior  to  field  experiments . 

When  studying  the  near  field  where  no  ionospheric  effects  are 
to  be  expected  and  where  the  earth  -  air  model  is  quite  sufficient 
for  simulating  the  field  conditions,  model  tests  are  suitable. 

Model  tests  in  conducting  media  thus  are  only  adequate  for  studying 
underground  antennas  (or  antennas  under  water)  and  are  certainly 
very  useful  as  long  as  the  used  model  factors  make  possible  the 
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simulation  of  the  antenna  dimensions  themselves.  Laboratory  tests 
with  antennas  whose  dimensions  are  approximately  equal  to  the  used 
wave  length  or  even  larger,  might  be  of  interest  in  the  sane 
way  as  investigations  of  the  directional  properties  of  antenna 
arrays  with  elements  having  distances  of  several  wave  lengths, 
but  they  cannot  be  applied  in  practice.  Their  results  therefore 
have  no  practical  importance  at  least  not  for  studying  VLF -waves. 

Summing  up  we  may  say  that  a  model  test  under  the  above  con¬ 
ditions  can  only  yield  useful  results  concerning  the  antenna 
properties  and  the  near  field.  For  reasons  of  better  understanding 

let  us  explain  the  expressions  near  field  and  far  field  used  in  the 

present  report.  In  order  to  decide  whether  a  point  of  observation 
is  situated  in  the  near  field  or  in  the  far  field  of  an  antenna, 
only  the  ratio  of  the  wave  length  to  the  distance  between  the 
site  of  observation  and  antenna  is  decisive  (in  contrast  to  the 
definition  of  some  authors,  the  antenna  size  is  not  important). 

Near  field  2tiR  «  X 

Far  field  2itR  »  A. 


4,2  Dipole  radiation  in  conducting  media 

The  field  of  an  electric  or  magnetic  dipole  embedded  in  a 
conducting  uniform  and  unbounded  medium  is  mathematically  described 
in  the  same  way  as  for  a  vacuum.  The  parameters  characterizing  the 
medium  are  to  be  replaced  by  the  corresponding  complex  values  in 
which  the  conductivity  is  taken  into  account.  Under  certain  circum¬ 
stances,  with  the  conductivity  being  sufficiently  high,  the 
displacement  current  may  be  neglected  and  simpler  expressions  may 
be  substituted  for  the  individual  parameters.  Thus,  for  example, 

+  ho  r'nmrvl  nn  rnnstmif  T^nrla 

~  w  ***x -  r“  '  r  -u  "  - -  - - -  —  .  --  .  — 
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iB*  =»  (i +  iwe))1/2  ^  (imps')'1'/2  =  (u>p<5/2)^'  2(l  +  i) 
=  B(l  +  i)  =  2n/\(l  +  i)  =  2n/\i  +  l/£  I 


(4,2,1) 


where  £  is  the  depth  of  penetration. 

The  wave  resistance  rj  =  E/2  ,  which  in  vacuo  is  (n/i)1/2,  also 
becomes  complex 


t-  £ 

O) 


(1  -  O  -  S<1  +  *>(4,2,2) 


Let  the  following  relations  be  valid: 


B  =  S€T 


A=  2(7r/fp<?)1/2  . 


f  ...  frequency 
6  • . .  conductivity 
A  » • •  wave  length 


=  2(n/f]it>y/*  .  (4,2,3) 

p  ...  permeability 
£  ...  dielectric  constant 
£  . . .  depth  of  penetration 


With  allowance  for  these  values,  the  following  expressions  are 
valid  for  an  electric  dipole  in  a  homogeneous  conducting  medium: 

F^=  (ll/4rc)  e“lBR/R(itovL  +  yj/R  +  1/Wrz  )  aiaS  (4,2,4a) 


=  (h/4ti)  e~iBR/R(iB  +  l/R)  sin 


(4,2,4b) 


For  the  magnetic  dipole  we  have: 


R/.-  (Kl/4n)e  iBR/R(iu>£  +  l/^R  +  l/iu>Rz  )  sin 
Ej»«  (Kl/47i)e"lBR/R(ifi  +  l/R)  sin  $- 


(4,2,5a) 

(4,2,5b) 


I  ...  current  intensity 
1  ...  antenna  length 


R,/£  ,y...  spatial  polar  coordihates 


magnetic  flux. 


Although  the  mathematical  derivation  of  the  expressions  (4,2,4) 
and  (4,2,5)  do  not  differ  from  those  for  a  dipole  vacuum,  their 
interpretation  and  further  requirements  resulting  from  them  are 
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i 

s 


quite  different.  First  of  ail  it  is  very  difficult  to  give  an 
adequate  definition  of  the  radiation  resistance.  In  air  the 
antenna  may  be  looked  upon  as  a  combination  of  the  ohmic 
resistance  and  the  radiation  resistance 


R  R0 

_ o  _ R 

-dm — i 


(A) 


The  radiation  resistance  is  obtained  by  calculating  the  ex¬ 
pression  for  the  total  energy  passing  through  any  spherical 
surface  that  envelopes  the  antenna,  and  dividing  it  by  the 
square  of  the  effective  input  current.  In  a  conducting  medium, 
the  energy  penetrating  a  spherical  surface  depends  on  the 
radius  of  the  sphere;  instead  of  (a)  we  now  have  to  consider 
the  antenna  and  its  neighborhoods 


RR 

CZ3 


(B) 


R  „  is  now  the  resistance  of  the  surrounding  medium.  This  means 
that  the  majority  of  the  energy  emitted  by  the  antenna  is  con¬ 
sumed  in  the  immediate  neighborhood  of  the  antenna  by  .the’ 
generation  of  Joulean  heat. 

Furthermore  it  can  be  shown  that  the  definition  of  the 
antenna  gain  used  in  vacuo  has  to  be  modified  for  a  con¬ 
ducting  medium.  In  the  far  field,  the  pattern  of  the  antenna 
field  lines  depends  on  the  choice  of  the  center  of  the  co¬ 
ordinate  system.  This  does  not  mean  that  the  structure  of  the 
field  depends  on  the  coordinate  system.  The  structure  itself 
is  invariant.  On  the  basis  of  the  exponential  attenuation,  the 
field  strength  at  a  point  of  observation  depends  on  the 
distance  from  an  antenna  point  in  a  different  manner  than  it 
does  for  a  vacuum. 

The  three  differences  between  dipole  radiation  in  vacuo 
and  in  a  conducting  medium  mentioned  herein  are  represented  in  M7"| 
Unfortunately ,  the  description  of  the  first  point,  i.e.  the 
calculation  of  the  radial  energy  flux,  shows  considerable 
errors  there.  An  exact  derivation  of  the  energy  phenomena,  a 
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comparison  of  the  electric  dipole  with  a  magnetic  dipole  which 
is  more  comprehensive  than  the  description  in  (jp)  as  well  as  a 
discussion  of  the  results  regarding  the  dependence  on  the  wave 
length  of  radiation  shall  be  given  in  what  follows. 

Prom  the  expressions  (4,2,4)  wc  obtain  the  following  ex¬ 
pression  for  the  radial  energy  flux: 


Pe  =  l/2Re  EjilJ= 


2v4rc‘ 

1/11x2  e 
2Mk; 


-r£R 
R5 
-2BR 


Hr  Rel>)>  +  b  +  fiB  +  - 


R 


/  0  2BS  2S 

(wp.B  +  +  ^vZ~  + 


R 


R 


(4,2,6a) 

(4,2,6b) 


i.e.  an  expression  which  depends  on  R  in  powers  higher  than  R  , 

-2£R 

which  also  contains  the  factor  e  .  In  the  case  of  vacuum  it 
reads : 


P®  =  1/2  ( 11/4*  )u>pfl/R* 


sm 


(4,2,6c) 


Integration  over  a  spherical  surface  of  radius  R  yields  the  to- 
tal  energy  W  propagating  through  this  surface. 


2 * 

\f  =  8n/3(ll/4*)|e-2BR/2j(wiifl  +  2SB/R  +  2S/R2  +  l/5R5). 

(4,2,6d) 


In  the  case  of  vacuum,  this  expression  is  reduced  to 


V/*  =  40k2  (llA)‘ 


(4, 2, fie) 


In  order  to  obtain  the  v/ave  length  dependence  of  (4,2,6d)  this 
expression  can  be  rewritten  by  means  of  (4,2,l)  and  (4,2,2): 

\f  =  (y)23e"47l//*R(4K2A5  +  2k//\25R  +  lA*R2  +  l/4  k*R5) 

=  (lL/A.)2je"^47l^R/3^(4K2//\  +  2k/r  +  A/r2  +  ^/4kR5). 


(4,2,6f) 
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where  no  term  is  independent  of  R.  For  R  »  A-  the  dependence 
on  R  is  reduced  to  the  exponential  factor,  yielding 

WR»A  =  4»V5(lVA)2c-"(4n//A)EM  . 

=  V/6air  e-(4"A)B/dA..  (4.2.68) 

Since  the  energy  propagating  through  a  spherical  surface  de¬ 
pends  on  the  radius  R,  the  radiation  resistance  cannot  be  de¬ 
fined  in  the  usual  way.  It  must  be  pointed  out  that  in  all 
expressions  (4.2.6)  and  in  the  calculation  that  follows,  passing 
to  the  limit  d  +  0  is  not  allowed,  since  the  condition  for 
neglecting  the  displacement  current  is  no  longer  fulfilled  for 
very  small  6,  6  may  increase  infinitely  in  the  same  way  as  A. 
Since,  however,  6  and  A  are  not  independent  of  each  other 
according  to  (4.2.3),  we  preferably  write 

6A  =  2(n/fvi6)1/2<$  =  2{nd/£4*lO~1)1/2  =  (10 1  6/f)1^2 
A  =  l/(lO"7f 6)1/2  . 


Thus,  (4.2,6g)  reads 


W. 


R»A 


=  W 


air 


^-'71(10"  7f6)1/2R 

(107  6/ f)1/2 


(4.2. 6h ) 


This  expression  holds  for  such  pairs  of  f  and  i  values  that  yield 
6/ 2nf  »£  .  The  energy  flux  in  the  magnetic  case  is  calculated 
from  the  expressions  (4.2.5)* 


Pm  =  -l/2  Re(Ej/>Hj)  = 

---  l/2(Kl/4Tt)2e''26R/R2  Re(i8  +  l/R)(iw£  +  1/qR  +  l/iu>pR2  j*sin%? 
=  l/2(Kl/4n)2e"23R/RZ^(B  +  l/R)sin2  J-  .  (4.2.7a) 


Here  too,  W  depends  on  R,  but  not  quite  as  strongly  as  in  the 


electric  case.  Integration  over  the  spherical  surface  of  radius 
R  yields: 

W®  =  (Kl/A  ) 2 1  e“ A/6i(  1  +  V4kR)  (4.2.7b) 

L- 

In  a  vacuum,  this  expression  is  reduced  to 

^  =  7i/3rl(Kl/>.)2 

Considering  the  limiting  case  R  »A.,  the  situation  is  quali¬ 
tatively  the  same  as  in  the  electric  case,  no  term  being  in¬ 
dependent  of  R;  quantitatively,  a  change  in  favor  of  the 
magnetic  dipole  takes  place  as  shall  be  shown  later  on. 

Yf£»A  .  (K1/a)W6  e-(4*/A)E  =  tf“alr  6OUoW/V^°-7«>1/2 

(4.2.7c) 

In  order  to  obtain  the  dependence  of  the  energy  fluxofcfand  f, 
we  write  the  following  expression  for  (4.2.6f)  and  (4.2.7b) 

We  "  ^air6"4^10  7i<r^  ^  /3 oUoV/f)1/2  £1  +  (l07/crf)1//2/2nR  + 

+  (lo7/r4/4tt2Ra  +llo7/f<r(io7/f£r)1//^/l6n5R5  (4.2.8a) 

L 

Wm  «  V/mair  60(l07<?/f)1/2  e"471^10  7fr^  ^  (1  +  (l07/cTf)1//2/27tR) 

(4.2.8b) 

Let  us  assume  that  there  exists  a  radius  Rq  so  that  the  energy 

radiation  of  both  dipoles  through  a  spherical  surface  having 

the  radius  R  is  equal.  Then  we  obtain  a  relation  between  the 
o 

energies  We.  and  VP.  and  thus  also  for  the  moments  II  and 
D  air  air 

K1  of  the  two  dipoles.  If  the  parenthetical  expressions  of 
(4.2.8a)  and  (4.2.8b)  are  denoted  with  h^,(R)  and  l^R),  re¬ 
spectively,  this  relation  reads  as  follows: 
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Under  this  assumption  we  consider  the  relation  if /if  at  a 


second  radius  R 


*7 V"  =  VJ®.r/l^.r|f/l800.107lJh1(R)/h2(R)  (4.3) 


Substituting  the  expressions  W  .  /ft  .  (4*2.9)  we  obtain 

•  air  air 


wVw*  =  h^R)  h2(Ro)/h2(R)h1(Ro)  =  (4.5b) 

±  -Li  “i  ++j}lPl  _I_  +  ip! L +  ip! 

;  2krJ  L1  Ui  2nR  df  2r2  tff  fdF  l6n3R3j 


i  +  n Z  _i_  1  H  ,  iliol  _i 

3  I  6f  2uR  _  L1  V  df  2nl 


1  107  1  107  H 107  _ l__2 

2*Rc  4*2R2  +  *f  16A5 

o  o 


(4.5c) 


For  R  » A.  the  following  expression  may  be  written: 


^ _ i)gf  ?hr0 _ 

W-  yio7  i  io7  i  iol  /I oZ  i 

1  *f  **0  4ti2R  2  *f  »e'f  16  A  5 

o  o 


llv  O' 


<  1  . 


(4.3d) 


The  expression  (4.3c)  increases  as  d  increases,  its  limiting  value 


being 


/We. 

lim  | 

\w“, 

air 


=  1  . 


This  shows  unambiguously  that  from  a  certain  d  =  .  onward 

mm 

and  at  a  given  frequency,  with  We.  /if1.  >  1,  the  magnetic  dipole 

air  air 

yields  better  results.  The  expression  (4.2.9)  is  larger  than 
unity  from  the  minimum  value  of  f  and  6  onward: 


Table  4.1 


f  (  g  t>o  ^  fmhno/ra^ 

min'  r 
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where  the  radius  Ro  has  been  put  equal  to  1  m.  Fig.  4..1  and  the 
corresponding  table  4*2  show  the  dependence  of  W6/??111  on  the 
frequency  in  the  range  f  =  1  -  100  kc/sec  for  three  different 
conductivities  ^  =  10  4,  10  ^  and  1  (mhos/m).  For  the  lowest 
curve,  the  ratio  We .  ftf1.  is  of  the  order  of  magnitude  10 
at  <1  =  10  mhos/m  it  is  10  and  at  cS  =  1  it  is  10  .  The  top¬ 
most  curve  along  which  V/6/^  lies  between  5*10~^  and  5  *10"" 4 
clearly  shows  that  even  though  the  energy  radiation  of  the 

electric  dipole  in  air  exceeds  that  of  the  magnetic  dipole 

5 

by  the  10  -fold,  in  a  conducting  medium  it  is  only  a  fraction 
of  this  value.  ' 


f(cps)  (^(mhos/m) 

10’ 

4-10' 

9.10^ 

16.10, 


36.10; 

64.10: 

81-10 

100.10 


3 


Table  4.2 
10~4 

0.79.10“^ 

-7-1  T 


3.1  .10 

7.1  .10 


1.2 

2.8 

3.4 

6.3 

7.8 


•  10 
•10 
>10 
-10 
-10 


-6 

-5 

-5 

-5 

-5 

-5 


-2 

10 

7.8.10“^ 
3.0.10**4 
6.9«10“i 
1.2*10“° 


2.7*10 

4.8*10 

5.9*10 

7.3*10 


-3 

-3 

-3 

-3 


5.9*10 
2.8-10 
5.8*10' 
9.5*10 
0.18 
0.28 
0.32 
0.30 


-3 

-2 

-2 

-2 


From  all  these  considerations  we  may  conclude  that  the 
magnetic  dipole  in  the  VLF  region  at  £  and  4  values  that  are 
of  practical  importance,  has  the  more  favorable  conditions 
for  the  generation  and  propagation  of  electromagnetic  fields 
in  a  conducting  medium. 

A  somewhat  different  approach  which  might  give  a  better 
explanation  of  the  physical  processes,  yields  the  same  result 
by  applying  the  theory  of  transmission  lines.  It  is  well  known 
that  the  dipole  field  may  be  looked  upon  as  being  a  combined 
derivative  of  a  scalar  potential  and  of  a  vector  potential  A, 
The  divergence  equation 

div  H  =  0 


for  the  electric  case  yields 


■f  -t 

H  =  curl  A  E  =  -imp  A  -  grad  V 
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i 


mohs/m 


001 
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In  the  magnetic  case,  ^  and  U  behave  like  current  and  voltage 

in  a  two-wire  cable  having  a  series  inductance  pH/m  and  a 

2 

shunt  inductance  of  H/m  and  a  shunt  capacitance  £.F/m. 


— -\-Vvv 


V\AA/ 


0» 


If  iu)fc*  =  iu)£  +6  is  substituted  for  iu>£  in  Eq.  (4.3.1),  then 
we  have 


§f  -  -(i»e  +^)v 


av 

3r 


-& 


lu>yi  + 


W 


2  S 

+  U?  £ 


HrT 


2£o) 


ICO 


— i — r 
to  t 


A  (4.3.3) 


for  the  electric  case.  The  equivalent-circuit  diagram  thus  is 
as  follows; 


VvA/v-i!- 


-Ww- 


H 


(E) 


It  contains  an  additional  series  resistance  of 
,  the  scries  capacity  thus  becomes 


26' 


(co262  +  J)t'  F 
2 to)2  m  f 

furthermore,  an  additional  shunt  resistance  ^  ^  is  present 
(in  all  cases,  is  to  have  the  amount  of  the  correspond¬ 

ing  physical  quantity  for  reasons  of  dimensioning,  because  A 
and  u  are  only  proportional  to  the  current). 

In  the  magnetic  case,  (4.3*2)  is  transformed  into 


4-13 


3 r 

dr 


-iwpU 


12 

3r 


2 

+  .  i 

iojy.r 


Equivalent  circuit  diagram: 


A/Wv 


•AA/Vv 


(4.3.4) 


(F) 


The  only  difference  is  an  additional  ohmic  loss  conductivity  <5" 

in  the  shunt.  The  comparison  of  (E)  and  (P)  shows  that  there 

exist  the  resistances  — -  and  - — r7~jr^~ — z — rr  without  any 

2  iwr  (6  +  u)£.  ) 

losses,  the  former  become  infinite  at  infinity  (as  expected), 

and  the  latter  disapoear  (at  infinity )$  furthermore  there 

2  6  ,  s 

exists  an  ohmic  resistance  ■  s""/  z ~g - rrr  in  (E)>  which  is 

r  (u>  e.  +  6  ) 

responsible  for  the  high  energy  loss  of  an  electric  dipole,  a.nd 
which  disappears  only  at  a  largo  distance  so  that  (e)  and 
(F)  there  become  identical.  This  is  the  physical  reason  for 
the  above  result,  namely  that  the  magnetic  dipole  is  more 
suitable  for  a  conducting  medium. 
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Summary 


In  the  present  report,  different  measuring  devices,  trans¬ 
mitting  and  receiving  anter.nas  (direction-finder  antennas)  for 
the  VLF  region  are  discussed.  The  measurements  of  propagation 
in  inhomogeneous  rock  are  mainly  concentrated  to  the  deter¬ 
mination  of  the  direction  of  the  field  vector  at  the  point  of 
reception.  The  problem  near  field  -  far  field  of  a  VLF  dipole 
is  also  discussed.  In  order  to  determine  the  electrical  pa¬ 
rameters  of  rock,  comprehensive  measurements  were  conducted 
in  solid  rock  and  in  rock  samples,  furthermore  the  frequency 
dependence  of  and  £.  was  studied.  The  influence  of  the 
gallery  configuration  on  the  result  of  a  £- measurement  was 
studied  in  detail.  The  theoretical  part  mainly  deals  with  the 
problems  of  excitation  of  a  VLF  wave  in  rock,  giving  a  detailed 
comparison  between  electrical  and  magnetic  dipoles. 
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